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1 Introduction 

Quantum statistics, as the name suggests, is the way that the quantum mechanical 
properties of particles influence their statistical mechanics. Two features of quantum 
mechanics are crucial for the difference between classical and quantum statistics: the 
indistinguishability of quantum particles and the 'finite room' that exists in a given 
domain of the phase space of a quantum system. 

The first feature can be incorporated 'by hand' into classical statistical mechanics 
(although there is no underlying principle for doing so) and this will resolve the 
Gibbs paradox, although by itself it would change little else. The second feature 
can likewise be incorporated in an ad hoc way, which will assign a unique value to 
the additive constant of the entropy of a classical system, but again without further 
physical impact. It is the combination of the two that produces nontrivial results, 
especially in the domain where quantum phenomena prevail, that is, high densities 
or low temperatures. 

The statistics of observed particles seem to be exhausted by the two well-known 
cases of Bose and Fermi. In this sense, the quest for more exotic possibilities seems 
to be a bit academic or even quixotic. It is, nevertheless, of interest to explore the 
possibilities as a way of obtaining a more lucid effective description of the dynamics 
of (otherwise normal) interacting particles. In this sense, we are trading one kind of 
intuition for another, hoping to make a profit in the bargain. 

There are several approaches to defining or implementing quantum statistics. 
They are roughly classified into 'fundamental' (group theoretic, operator etc.) and 
'phenomenological' (statistical). While exploring them we must keep in mind that, 
at the end of the day, we should provide a realization, that is, a system which 
is (more or less) physical and exhibits these statistics through its dynamics. The 
less the number of space dimensions, the more easy it is to find such realizations. 
The concept itself of statistics, actually, depends crucially on the dimensionality 
of space. Topological considerations allow for generalized abelian statistics in 2+1 
dimensions [|^, ^ ^ Q. Such particles, called anyons, have found use in the physics of 
the quantum Hall effect and (probaby) high-temperature superconductivity. There 
are excellent reviews on anyons |^ and we will not deal with them here. 

In these lectures the bulk of the material will concern one- dimensional situations, 
although some considerations in principle applicable to any dimension will also be 
included. The purpose is to present enough of the concepts and material on the 
subject to spur the interest of the reader for further study. You are warned that 
all formulae in these lectures are simply mnemonics of the correct results: mistakes 
abound! (The written version is only marginally better.) One standard exercise 
proposed here is to rederive all results, fixing the signs and factors in the process. 
Hopefully substance survived the sloppyness. 

So, right to the task. 



2 Permutation group approach 

As a way of introduction and covering standard ground, we will study the context in 
which statistics initially arose, namely the permutation group in quantum mechanics 
1^-0]. The main emphasis will be in connecting to the path integral and pointing 
the possibilities for generalizations within this approach. 
The main themes of this lecture are: 

a) Eliminating the 'gauge' multiplicity of states originating from particle permuta- 
tions. This will, in general, lead to the appearance of internal degrees of freedom 
that are the residue of 'identities' (distinguishability) of particles. 

b) Implementing the path integral for identical particles. Their indistinguishability 
calls for including paths that lead to any permutation of the final position. Deter- 
mining the right weights of these sectors (which are not necessarily phases) will fix 
the statistics. 



2.1 Realization of the reduced Hilbert space 

We will work in the first-quantized picture, in which the number of particles is 
fixed but the particles are nevertheless taken to be indistinguishable. This is to 
be contrasted to the weaker statement that the particles are merely identical. To 
fix the ideas, let us call {xj} the coordinates of these particles, with i = 1 . . . N 
enumerating the particles, and assume, for now, that there are no internal degrees 
of freedom. The Xi can be in a space of any dimension. The Hilbert space of this 
system is spanned by the position eigenstates \xi, . . .xn >= \x >. Let us also define 
the operators P E Sjy that permute the quantum numbers of the particles. There 
are A^! such operators forming the permutation (or symmetric) group of A^ objects 
Sn. The collection of the above states for a set of distinct Xj transforms in the 
A^!-dimensional regular (defining) representation of S^ 

P|X >= |PX >= |xp(i),...Xp(Af) > (1) 

Notice that if any of the coordinates Xi coincide the above is not the full defining 
representation any more. The set of such states, however, is of measure zero (the 
coordinate space assumed continuous) and thus they can be safely ignored. (We 
assume that there are no interactions singular at coincidence points that might 
dynamically make such states of nonzero weight.) 

If the particles are identical then the hamiltonian treats them all equally and 
thus we have 

[H, P] = (2) 

for any P. If however, the particles are also indistinguishable, this means that there 
is no physical way to ever be able to tell the particles from each other. Thus, not 
only the hamiltonian, but all conceivable perturbations of this hamiltonian must 
commute with the permutation operators. (These perturbations would correspond 



to various devices, dies and other tricks that we would attach to the system in our 
effort to tell the particles from each other.) In short, the permutation operators 
must commute with all physical observables A: 

[A, P] = (3) 

When such an operator P commutes with all (physical) operators in a system we 
have a superselection rule: there is no possible transition or interference between 
states corresponding to different eigenvalues of this operator and thus we can choose 
a subspace with a fixed eigenvalue and stay with it for ever. (This is also called 
'picking a theta- vacuum' in particle gauge theories.) In our case, there is a whole set 
of operators commuting with everything else (the group Sn) and thus the generahzed 
'theta- vacuum' states consist of irreducible representations (irreps) of Sn- Further, 
states that are connected to each other through the action of operators P cannot 
be distinguished by any physical operator and therefore must be identified as a 
unique physical state. In other words, the permutation group must be viewed as 
a discrete gauge group producing unphysical transformations, that is, copies of the 
same physical system. We must, therefore: 

• Project on a subspace corresponding to a particular irrep R of Sn 

• Identify elements within this subspace grouping in the same irrep as a unique 
physical state. 

The above procedure and choice of R corresponds to a choice of statistics. 

We will implement this in the coordinate representation \xi,...xn >= \x >. 
Projecting the Hilbert space to an irrep R of Sn amounts to keeping only linear 
combinations of states within this multiplet that transform in R, that is, 

\a;x>='^Ca{P)P\x> , a=l,...dji, dR = dim{R). (4) 

p 

where the sum is over all elements of the permutation group and Ca{P) are appro- 
priately chosen coefficients. If we denote with Rab{P) the matrix elements of the 
permutation P in the representations R, 

P\a,x>=Y.RabiP-')\b,x> (5) 

b 

(The appearance of P~^ above is necessary so that successive transformations com- 
bine in the right order.) 

The defining representation decomposes into irreducible components, classified 
by Young tableaux, each appearing with a certain multiplicity. Should we keep only 
one irrep out of each multiplicity or the whole multiplet? To decide it, note that if 
instead of the base state \x > for the construction of the states \a,x > we choose 
a different permutation Po\x >, then although the new states \a, PqX > constructed 
through (^) still transform in the irrep R, in general they are not linear combinations 
of \a,x > but rather span a different copy of R. Since we can continuously move 



in the configuration space from \x > to Po\x >, we conclude tfiat we must keep 
all irreps R witliin eacli multiplet. (In otlier words, altfiougfi for eacli point in the 
Hilbert space \x > this multiplet is reducible, the fiber of these representations over 
the Hilbert space is connected and irreducible.) 

To realize explicitly the above, we construct the states 




\ab,x>=d^y^R,,iP)P\x> (6) 



Using the group property of the representation R{Pi)R{P2) = R{PiP2), we deduce 
that under the action of the group Sn and under change of base point x the above 
states transform as: 

P\ab, X >= J2 RaciP'^)\cb, X > , \ab, Px >= J2 Rcb{P'^)\ac, X > (7) 

c c 

Thus we see that the first index in these states labels the different elements of a 
single irrep R, while the second index labels the different equivalent irreps in the 
multiplet. Since both indices take cIr values, we recover the standard result that 
each irrep of Sn is embedded in the defining representation a number of times equal 
to its dimension. 

Consider now the matrix element < ab,x\A\cd,y >, where A is any physical 
operator, that is, any operator commuting with all elements P of S^- Substituting 
the definition (^) and using the unitarity of P {P'^ = P~^) and of R {R*^^{P) = 
Rba{P~^)) we obtain, after a change in summation variable, 

< ab,x\A\cd,y>= |^ ^ R,,{P')RUP-')Rcd{P) < x\AP'\y > (8) 

^^ ■ P,P',e 

Using further the orthogonality (Shur's) relation (see, e.g., |15 



J^Rab{P)Rcd{P-') = ^<5aAc (9) 

p UR 

we finally obtain 

< ab, x\A\cd, y>=J2 ^acRbd{P) < x\A\Py > (10) 

p 

It is clear from (^) that there is no possible transition between states with different 
first index. Further, for states with the same index, the amplitude is independent 
of the index. Thus, the first index a in the state \ab,x > propagates trivially. Since 
this is the index that corresponds to the different but physically equivalent states 
within each irrep R, we conclude that the required projection of the Hilbert space 
to the physical subspace amounts to simply omitting this index from all states. 
(That is, freeze this index to the same fixed value for all states of the theory; no 



physical quantity will ever depend on the choice of this value.) On the other hand, 
the second index, corresponding to different equivalent irreps, does not propagate 
trivially and must, as argued before, be kept. We are led therefore to the physical 
states \ba,x >— *> \a,x >. 

Let us first choose A = 1. Then ( p!0D provides the overlap between the states 

< ab, x\cd, y>=Y. SacRbd{P)S{x - Py) (11) 

p 

For X in the neighborhood oiyitisP = l which contributes to the normalization, for 
which Rbdi^) = ^bd and we recover the standard continuous normalization between 
the states. 

Choose, now, A = e~*^*, where H is the hamiltonian, and thus find the propa- 
gator G{ab,x;cd,y\t) between the states of the system. Projecting to the physical 
states \ba,x >— *> \a,x > we obtain the propagator 

Guia, X- 6, y\t) = Y. R^,{P) G{x, Py; t) (12) 

p 

where G{x,Py;t) = < x\e~^^^P\y > is the usual many-body propagator. We 
note that, due to the transformation property (0), the states \a,Px > are linear 
combinations of states \a,x >. Therefore, projecting down to the physical subspace 
corresponding to R amounts to trading the original A^! copies of physically equivalent 
states \Px > for a number dn of global internal degrees of freedom for the system, 
labeled by the index a. 

2.2 Path integral and generalized statistics 

It is now easy to write down the path integral corresponding to identical particles 
quantized in the i?-irrep of S^. G{x, Py; t) can be expressed as an A^-body path 
integral in the standard way, with particles starting from positions Xi and ending 
in positions Pyi = yp{i). Since all permutations of particle positions are physically 
equivalent, (12) instructs us to sum over all sectors where particles end up in such 



permuted positions, weighted with the factors Rab{P) depending on the internal 
degrees of freedom of the initial and final states. From (^, |ll]) we can write the 
completeness relation within the physical subspace 

f d^x 
Ir = J —^Y\a,x><a,x\ (13) 

and with the use of (^) it is easy to prove that the above path integral is unitary, 

that is, 

r d^v 
y ^ E G{a, x; b, y\t) G{b, y; c, z\t') = G{a, x; c, z\t + t') (14) 

We can extend this construction to a more general class of statistics, which 
includes the cases of the so-called 'parabosons' and 'parafermions' as special cases. 
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The generalization consists in allowing more than one irrep to appear in the Hilbert 
space, and keeping more than one state in each irrep. This seems unmotivated 
and against the spirit of the reduction by the full superselection rule, but it is 
consistent. It could mean, for instance, that the particles have some hidden internal 
degrees of freedom accounting for the extra degeneracy, which are invisible to the 
present hamiltonian but may become dynamically relevant later. (This, then, would 
augment the space of physical operators in the theory.) The most general situation 
is that we include Cr states from each irrep R. The statistics is fixed by the set of 
non-negative integers Cr. 

The generalization of the results to this more general situation is straightforward. 
The internal degree of freedom now takes values A = {R, a, a), where a = 1, . . .dR 
labels the inequivalent /2-irreps and and a = 1, . . . Cr labels the states kept within 
each irrep. (a, a) constitute the internal degrees of freedom within each irrep. So, 
overall, A takes J2r CRdR different values. The propagator (and corresponding path 
integral) is obviously 

Gs{A, X- B, y\t) = Y. S{P)ab ^(x, Py; t) (15) 

p 

where 

SiP)AB = Sr^,r^ 5^p {RA)ab{P) (16) 

The definition of parabosons (parafermions) of order p is that we include once all 
irreps with up to p rows (columns) in their Young tableaux. We have, therefore, 
Cr = 1 for such irreps and Cr = for the rest. We note that the irreps for 
parafermions are the duals of those for parabosons (the dual of a tableau is the 
tableau with rows and columns interchanged). In an appropriate basis, the repre- 
sentation matrices of dual irreps R, R are real and satisfy 

RabiP) = i-lfRabiP) (17) 

where (—1)^ is the parity of the permutation. We arrive then at the relation between 
the weights for parabosons and parafermions of order p: 

Spf{P)ab = {-lfSpB{P)AB (18) 



This extends a similar relation for ordinary fermions and bosons |jT^, for which there 
are no internal degrees of freedom and Sb{P) = 1- 

From the path integral we can evaluate the partition function, by simply shifting 
to the euclidean periodic propagator Ge{P) = e~^^ and summing over all initial 
and final states, with the measure implied by ([T3|). Given that 



J2RaaiP)=tTRiP)=XRiP) (19) 



we get the expression in terms of the characters of Sn 



Zs{T) =J^J:S{P) < x\GE{f3)\Px > (20) 



where 

S{P) = Y.CrXr{P) (21) 

R 

We also note the inversion formula 

Cr = ^T.S{P)xr{P) (22) 

which allows to recover Cr if we know S{P) for all permutations. The interpretation 
in terms of a periodic euclidean path integral is obvious. The characters Xr{P) ^-'^s 
a set of integers, and thus the "statistical factors" S{P) weighing each topological 
sector of the path integral are (positive or negative) integers. In the case of para- 
bosons of any order p, however, we note that the statistical weights are positive (or 
zero) integers. The ones for parafermions can be either positive or negative, as given 
by 

SpF{P) = i-lfSpBiP) > SpBiP) > (23) 

We do not have a general formula for SpB{P) for arbitrary p. For non- interacting 
particles the partition function can be expressed in terms of the characters of SU{N) 
(Shur's functions) |T^, [1^ . 



We note here that the case of distinguishable particles ("infinite statistics" 10, 



T8| , 0), is also included in this formalism, by accepting all states in each irrep, that 
is, Cr = cIr. Since R appears exactly (Ir times in the defining representation of Sn, 
S{P) above becomes the trace of P in the defining representation. But all P 7^ 1 are 
off-diagonal in the defining representation, so we get Sinf{P) = N\ 6p^i, recovering 

Z{T)= fd^x<x\GE{(3)\x> (24) 

for distinguishable particles. 

2.3 Cluster decomposition and factorizability 

Parastatistics particle obey the cluster decomposition principle, in the sense that the 
density matrix obtained by tracing over a subset of particles which decouple from 
the system can be constructed as a possible density matrix of the reduced system of 



remaining particles [g| |T2|. From (^^, however, we see that the partition function 
of two dynamically isolated sets of particles A^^i and N2 does not factorize into the 
product of the two partition functions, since the statistical weights S{P) in general 
do not factorize into S{Pi)S{P2) when P is the product of two commuting elements 
Pi and P2. Equivalently, this means that the occupation degeneracy D{pi,p2, . . .) 
oi pi,p2, . . . particles occupying a set of levels €1,62, ■■ ■ does not factorize into the 
product of individual occupation degeneracies for each level e^. There is an effective 
'coupling' between the particles. This has important physical implications. If the 
two sets of particles are totally isolated, it does not make sense to evaluate the 
partition function of the total system, since the statistical distribution can never 



relax to the one predicted by that partition function. The individual partition 
functions of the subsystems are the relevant ones. If, however, the two sets are 
only weakly coupled, then initially each set will distribute according to its reduced 
partition function, but after some relaxation time (depending on the strength of 
the coupling between the two sets) they will relax to the joint distribution function, 
which, we stress, will not even approximately equal the product of the individual 
ones. Thus, cluster decomposition holds in an absolute sense but fails in a more 
realistic sense. In contrast, fermions and bosons respect cluster decomposition in 
both senses. 

We summarize by pointing out that the most general statistics of the type ex- 
amined here is parametrized by any of three possible sets of numbers. The first is, 
as just stated, the number of states Cr accepted for each irrep Roi Sjy. The second 
is the statistical weights S{P) appearing in the partition function (euclidean path 
integral). These weights are invariant under conjugation oi P ^ QPQ~^. Finally, 
we could use the degeneracy of a many-body occupancy state, call it D[pi], as our 
definition. (This was the starting point of the first known attempt to generalized 
statistics, the 'Gentile statistics.' I^O]) These three sets contain as many elements 
as the partitions of A^ and are all equivalent. 

What are the restrictions or criteria to be imposed on the above parameters? 
The first one is unitarity, that is, the existence of a well-defined Hilbert space with 
positive metric. This requires that C/j be non-negative (no negative norm states) 
integers (no "fractional dimension" states). The other will be what we call "strong 
cluster decomposition principle" or "factorizability," that the partition function of 
isolated systems factorizes. This is a physical criterion, rather than a consistency 
requirement. To summarize: 

• Unitarity: Cr non-negative integers 

• Factorizability: S{P) = Iln&cycies{P)Sin) or D[pi\ = UiD{p,) 
Factorizability, in particular, implies that the grand partition function in the case of 
noninteracting particles will be obtained by exponentiating the sum of all connected 
path integrals [P a cyclic permutation of degree n) with weights S{n)/n {1/n is 
the symmetry factor of this path integral, corresponding to cyclic relabelings of the 
particles). The grand partition function will factorize into a product of partition 
functions for each level e^. Thus, S{n) are cluster coefficients connected to D{n) in 
the standard way 

f:D{p)z^ = exp(f:^zA (25) 

p=0 \n=l ^ / 

The above formula, in fact, provides the easiest way to relate D[pi\ and S{P) in the 
general case (no strong cluster property) but we will not enter into this. 

If we assume that S{1) = D{1) = 1, then it is easy to verify that the only solution 
of the above two criteria is ordinary fermions and bosons. The situation is different, 
however, when S'(l) = -D(l) = g > 1 (this would mean, e.g., that the particles come 
a priori in g different "flavors"). The possibilities are manifold. All these generalized 
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statistics share the following generic features: 

• The degeneracy of the state where n particles occupy different levels is g". 
(Indeed, £>(!, 1, ... 1) = D{lY = g".) 

• If state A can be obtained from state B by 'lumping' together particles that 
previously occupied different levels, then D{A) < D{B). (E.g., -D(3) < D{2,1) < 
D(l,l,l).) 

The above possibilities include the obvious special cases of gi bosonic flavors and 
g2 fermionic ones (gi + g2 = g), for which S{n) = gi — (— l)"g2, along with many 
other. As an example, we give the first few degeneracies for many-particle level 
occupation for all statistics with g = 2: 

D(l) = 2, D{2)=4, D(3) = 8 

D(l) = 2, 1^(2) =3, D(3) = 6,5,4(5 + 5) 

D{1) = 2, D{2)=2, D(3) = 4,3,2(5 + F),l,0 

D(l) = 2, D{2) = 1, D(3) = 2,1,0(F + F) 

D{1) = 2, D{2)=0, D(3) = 

The specific choices denoted by B + B, B + F and F + F are the ones corre- 
sponding to two bosonic, one bosonic and one fermionic, and two fermionic flavors 
respectively. The topmost statistics could be termed "superbosons" and the bottom 
one "superfermions" of order 2. We also remark here that the "(p, g)-statistics" 
introduced in |jl2| can be realized as particles with p bosonic and g fermionic flavors, 
where we identify each multiplet transforming irreducibly under the supergroup 
SU{p, g) as a unique physical state. 



Finally, we comment on 'Gentile statistics' pO|. The rule is simply that up to p 
particles can be put in each single-particle level. This corresponds to D{n) = 1 for 
n < p, and D{n) = otherwise. This has been criticized |0 on the grounds that fix- 
ing the allowed occupations for each single-particle state is not a statement invariant 
under change of single-particle basis. It should be clear from this lecture that any 
statistics satisfying the unitarity requirement is consistent and basis- independent. 
Therefore, Gentile statistics must violate unitarity. Indeed, it is easy to check that 
all weights C/j for such statistics are integers (this is generic for all statistics with 
integer D{n)), but not necessarily positive. In the specific case of p = 2, e.g., where 
up to double occupancy of each level is allowed, the degeneracies of each irrep of Sn 
(parametrized, as usual, by the length of Young tableau rows) up to A^ = 5 are 

C, = C21 = C22 = C221 = 1, Cm = Cnii = C2111 = -1, else Cr = (26) 

We see that representations 111, 1111, 2111 correspond to ghost (negative norm) 
states and their effect is to subtract (rather than add) degrees of freedom. 

In conclusion, we see that the permutation group approach to generalizing statis- 
tics gives results of rather limited interest. The resulting statistics can be though 
of as particles with internal degrees of freedom, with perhaps some superselection 
rule that forbids the manifestation of all possible internal states. To find something 
more exotic we must consider other approaches. 
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3 One-dimensional systems: Calogero model 

Let us now specialize in one dimension and see in what different angle we can 
approach the problem of statistics. The peculiarity of a lineal world is that, kine- 
matically, particles cannot exchange positions without "bumping" onto each other. 
This has tangible physical consequences and makes the notion of statistics in one 
dimension rather special. Let's see: 

A. Since configurations related by particle permutation are 'gauge' copies of 
each other, we could 'gauge fix' and restrict ourselves to only one gauge copy. In 
one dimension the configuration space of A^ indistinguishable particles breaks into 
N\ sectors classified by the ordering of the coordinates: Xi < ... < xn and its 
permutations. Restricting to one sector we are faced with the problem of boundary 
conditions on the wavefunction at the boundary of this space, when two or more 
coordinates become equal, so that we preserve hermiticity. The choice of boundary 
conditions can be interpreted as a choice of statistics. We mention two possibilities: 

a) A linear local boundary condition: 

V' + Aa„^|,„,,,,,^ = (27) 

where 9„ is the normal derivative at the boundary |^. Clearly A = (Neumann) 
corresponds to fermions; we can analytically extend ip in the other sectors in an 
antisymmetric way. Similarly, A = oo (Dirichlet) corresponds to bosons. Any other 
choice would be some intermediate statistics. A, however, introduces a length scale 
and is not a very satisfactory statistics definition. At any rate, this system is equiv- 
alent to bosons with a delta-function two-body interaction of strength proportional 
to A^i. 

b) Fix the analytical behavior of ip at the boundary as 

ij r^x^ as x ^ (28) 

This way the probability current would scale as ipdnip ~ £a;^^~^ and would vanish if 
£ > I or f = 0. This would define the statistics through the dimensionless parameter 
(,. Notice that this behavior would require a two-body potential behaving like ' ^ ' 
near coincidence points. This is the first ghmpse at the inverse-square potential 
arising in the statistics context. 

B. We could, alternatively, examine the scattering phase between particles. At 
high relative energies the two-body scattering phase approaches the values 

^ ^ for bosons, 6 ^ tt for fermions (29) 

irrespective of the details of their interaction (provided it is not too singular). The 
corresponding many-body phase is the sum of two-body ones. So, if we see a system 
where this phase goes to the value in we can interpret it as one with generalized 



statistics of order i 21 . 
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As we shall see, this definition also leads to the inverse-square potential. There 



are hints from algebraic considerations also pointing to this type of potential |22, 



531, E3 . This is enough motivation, at any rate, to examine particles with that type 



of two-body interaction potentials. This is known as the Calogero model, or, more 



completely, as the Calogero-Sutherland-Moser (CSM) model p5|, ^ ^ |. 



3.1 The Calogero-Sutherland-Moser model 

This is a system of nonrelativistic identical particles on the line with pairwise inverse- 
square interactions. The basic hamiltonian is 

This is the 'free' (scattering) Calogero system on the line [^. The particle masses 
m have been scaled to unity. This is the only scale-free two-body potential that 
one can have (the potential, quantum mechanically, scales like the kinetic term). 
An external harmonic oscillator potential can also be added to confine the system 
without spoiling its features. This is the harmonic Calogero model. Alternatively, 
one could consider a periodic version of the system. The particles now interact 
through infinitely many periodic copies of themselves and the two-body potential 
becomes 

^(^)= £ f^ , L^^2 = T^^ (31) 



(x + 27rn)2 ^2 



sm| 



(we scaled the length of the box to 27r). This is the Sutherland model |26|. There 



are other versions of this class of models that will not concern us here. The classic 



report |28] covers these systems in detail. See also [29| for recent mathematical 
work. 

Classically the coupling constant g should be positive to ensure particles are not 
'sucked' into each other. Quantum mechanically the uncertainty principle works in 
our favor and the minimum allowed value ior g is g = —^ (we put henceforth ^ = 1). 
For later convenience, it is useful to parametrize g in the fashion 

g = i{i- 1) (32) 

in which case the minimum value is naturally obtained for i = ^■ 

The above system is integrable, which means that there are N integrals of motion 
in convolution, that is, A^ functions on phase space with vanishing Poisson brackets: 

{In,Im} = 0, n,m = l...N (33) 

For the scattering system (no external potential) Ji is the total momentum, I2 is the 
total energy amd the higher /„ are higher polynomials in the momenta also involving 
the two-body potentials. We will not belabor here their form nor prove integrability 
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at this point, since this is beyond our scope and will, at any rate, be dealt with in 
later lectures. We will only state without proof the qualitative features relevant to 
our purposes. 

The key interesting property of the above model, that sets it apart from other 
merely integrable models, is that, both classically and quantum mechanically, it 
mimics as closely as possible a system of free particles. Let us first look at its classical 
behavior. The motion is a scattering event. Asymptotically, at times t = ±00, the 
particles are far away, the potentials drop off to zero and motion is free. When they 
come together, of course, they interact and steer away from their straight paths. 
Interestingly, however, when they are done interacting, they resume their previous 
free paths as if nothing happened. Not only are their asymptotic momenta the same 
as before scattering, but also the asymptotic positions (scattering parameters) are 
the same. There is no time delay of the particles at the scattering region. The only 
effect is an overall reshuffling of the particles. Thus, if one cannot tell particles from 
each other, and if one only looks at scattering properties, the system looks free! 

This behavior carries over to quantum mechanics. The asymptotic scattering 
momenta are the same before and after scattering. The fact that there is, further, 
no time delay translates into the fact that the scattering phase shift is independent 
of the momenta. Thus it can only be a function of the coupling constant and the 
total number of particles. It is, actually, a very suggestive function: 

e.^.^[^i. (34) 

Thus the phase is simply in times the total number of two-body exchanges that 
would occur in the scattering of free particles. Clearly the case £ = would corre- 
spond to free bosons and i = 1 to free fermions (for these two values the potential 
vanishes and the system is, indeed, free). For any other value we can interpret this 
system as free particles with generalized statistics £. 

A word on the permutation properties of this system is in order. The inverse- 
square potential is quantum mechanically impenetrable, and thus the 'ordinary' 
statistics of the particles (symmetry of the wavefunction) is irrelevant: if the particles 
are in one of the A^! ordering sectors they will stay there for ever. The wavefunction 
could be extended to the other sectors in a symmetric, antisymmetric or any other 
way, but this is irrelevant for physics. No interference between the sectors will ever 
take place. All states have a trivial A^! degeneracy. (That is to say all states in all 
irreps of S^ have the same physical properties.) Permutation statistics are therefore 
irrelevant and we can safely talk about the effective statistics as produced by their 
coupling constant i. 

The behavior of the wavefunction near coincidence points in as in (^). This 
system, thus, satisfies both our boundary condition and scattering phase criteria for 
generalized statistics. 

Let us also review the properties of the confined systems. In the presence of an 
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external harmonic potential of the form 



V = Y. \^'^\ (35) 

i 

the energy spectrum of a system of uncoupled particles would be 

E = —u + Y^ UiU) (36) 

^ i 

The Hi are nonnegative integers satisfying 

Til < ... < riN for bosons (37) 

rii < . . . < rijq for fermions (38) 

The actual spectrum of this model is 

E = -tu + £^— '-u + Y^Tiiu; (39) 

with TLi being 'pseudo-excitation numbers' obeying bosonic selection rules: rij < rij+i. 
Defining the 'quasi-excitation numbers' 

n^ = ni^(i- \)i (40) 

we can check that the expression of the spectrum in terms of the fii is identical to 
the free one ( pG] ) but with the quantum numbers now obeying the selection rule 

fii < Ui+i - i (41) 

This is a sort of exclusion principle that requires the particle quantum numbers to 
be at least a distance £ apart (as contrasted to 1 for fermions and for bosons). 
Again, a generalized statistics interpretation is manifest ||3CI|| . 

Let us clarify that the above numbers rij are no more integers. They do, however, 
increase in integer increments. The rule is that the ground state is determined by 
the minimal allowed nonnegative values for n, obeying (^I]) while the excited states 
are obtained by all integer increments of these values that still obey (^Tj). 

The periodic (Sutherland) model has similar properties. Its spectrum is 

E = T. \k^ + ^E(% - kd + a' ^^^l'^^ (42) 

with the 'pseudo-momenta' hi being integers satisfying bosonic rules: ki < /cj+i. This 
looks rather different than the corresponding free expression (for ^ = 0). Defining, 
however, again 'quasi- moment a' 

P^ = h + i(^^-^^^^ (43) 
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we can check that the expression for the spectrum becomes 

E = Y.\p" (44) 

that is, the free expression. The quasi-momenta satisfy 

Pi < Pi+i - ^ (45) 

that is, the same selection rule as the n^ before. Again, we observe a generalization of 
the fermionic and bosonic selection rules corresponding to statistics I. The ground 
state is the minimal (nearest to zero) numbers satisfying ( ^5]) while excitations 



correspond to integer increments (or decrements) of the pi still satisfying (^5]). 

Let us note that the above rule for i = \ reproduces the fermionic spectrum of 
particles with periodic boundary conditions for odd A^ and anti-periodic ones for 
even iV: in the odd (even) iV case the momenta are quantized to (half-) integers. 
This has a natural interpretation: when we take a particle around the circle it goes 
over A^ — 1 other particles. If we require the phase shift of the wavefunction in this 
process to agree with the minus signs picked up from the A^ — 1 fermion exchanges we 
recover the previous rule. We stress that, for free particles, this is not a consistency 
requirement but rather an aesthetic rule. At any rate, this is what the Sutherland 
model chooses to do! 

In conclusion we see that the CSM model can be though of as a system of particles 
obeying generalized statistics. This manifests in terms of the boundary conditions of 
the wavefunction, the scattering phases and, most significantly, through a peculiar 
'level repulsion' of their quantum numbers generalizing the Fermi exclusion principle. 

3.2 Large-N properties of the CSM model and duality 

Let us examine, now, the properties of the CSM model as the number of particles 
grows large. At zero temperature, a non-interacting fermion system would form 
a Fermi sea. The corresponding 'Fermi surface' in one dimension degenerates to 
points. For the system in an external harmonic potential there is just one point 
corresponding to the highest excitation np = N — \. For the free periodic system 
we would have two Fermi momenta at ipi? = ±^^^--!-. Excitations over this ground 
state are, then, conveniently classified in terms of particles (a filled Fermi sea with 
an isolated particle above or below) and holes (a filled sea with one unoccupied state 
inside it). 

Interestingly, the CSM model presents a similar picture. The qualitative features 
of both the Calogero and the Sutherland model are similar, so we pick the latter as 
most closely representing a gas of free particles in a box. From (^5]) above we see 
that the ground state also forms a 'quasi-Fermi sea' (or should we call it a 'Luttinger 
sea'?) with Fermi levels rescaled by a factor t. pp = f-^^-^- Its minimal excitations 
are analogous to the ones of a Fermi sea, but not quite: 
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• A particle would be an isolated occupied quasimomentum above or below a 
completely filled sea of quasimomenta. Particles 

-are excited in units of 1 (the increments of their quasimomentum). 

-take up a space C, in quasimomentum (since they cannot be 'packed' closer than £ 

units apart). 

• A hole would be an isolated empty space inside an otherwise occupied sea. 
Interestingly, the minimal such excitation is not obtained by removing one particle 
from the sea, but rather by incrementing all quasimomenta of the sea above the 
place where we want to create the hole by one unit. Holes 

-are excited in units of L Indeed, since the distance of quasimomenta in the sea is 
£, the possible positions of the hole are at distances C, apart. 

-take up a unit space in quasimomentum. Indeed, incrementing all quasimomenta 
above a given place in the sea by two units creates two holes in that place, and so 
on; by locally reshuffling quasimomenta we can then separate these holes. 

Note that holes are not antiparticles. Removing a particle for the sea creates a 
gap of i spaces and, from above, C, holes. The correspondence is 

1 particle ~ — £ holes (46) 

We already observe a sort of duality between the two types of excitations. This can 
be summarized as 

particle <-> hole , £<->-, p ^^ ip (47) 

Under the above, the spectrum of excitations of the model remains invariant. This 
is the simplest manifestation of a coupling-constant duality that goes over to the 



correlation functions and Green's functions of the model ||31], ^ |33|, 0]. Obviously, 
this duality is spoiled by nonperturbative effects, since holes are confined within the 
sea while there is no 'ceiling' for particle excitations. 

This concludes our brief 'physicist's tour' of the CSM model. We are now going 
to examine the most useful ways of analyzing this system. 



4 One-dimensional systems: Matrix model 

Identical particles can be formulated in terms of their phase space variables modulo 
permutations. This leads to the symmetric group approach of section 1. We could 
try other approaches. For instance, we could formulate them in terms of the eigen- 
values of an A^ X A^ matrix. There is no a priori ordering of these eigenvalues, so this 
certainly encodes identical particles. It is clear that the permutation symmetry of 
the problem has been promoted to the continuous symmetry of unitary conjugations 
of this matrix, which leaves the eigenvalues intact. This should be viewed as a gauge 
symmetry and would open the road to yet another definition of statistics. 
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4.1 Hermitian matrix model 

With this introductory remark as motivation, let us examine a matrix model that 
parallels as closely as possible particle mechanics. The kinematical variable is a 
hermitian N x N matrix M and the lagrangian reads 

C = tTi-M^-V{M)\ (48) 

V{x) is a scalar potential evaluated for the matrix variable M. 

Clearly the above has a time-translation invariance which leads to the conserved 
energy 

H = tTl^M^ + V{M)\ (49) 

Moreover, the action is invariant under time-independent unitary conjugations of 
the matrix M: 

M -^ UMU-^ (50) 

This nonabelian SU{N) symmetry leads to the conserved hermitian traceless matrix 

J = i[M,M] (51) 

where [ , ] denotes ordinary matrix commutator. These are the 'gauge charges' that, 
when fixed, will determine the 'statistics' of the model. But let us first examine the 
system classically. We are interested in the dynamics of the eigenvalue of M, so we 
parametrize it as 

M = UAU-^ (52) 

where U{t) is the unitary 'angular' part of the matrix and A(t) = diag{xi, . . .Xn} 
are the eigenvalues. Clearly the conserved quantity J has to do with invariance 
under 'rotations' of the angular part of M and thus corresponds to the 'angular 
momentum' of U{t). We define the 'gauge potential' 

A = -U-^U (53) 

M, J and the lagrangian C become, in this parametrization, 

M = f/ (A + [A, A]) f/-^ (54) 

J = iU{[A,[A,A]])U-^ = UKU-^ (55) 

C = tr|iA2 + [A,A]2-\/(A)| (56) 

1 TV 1 TV 

1=1 *ii=i 



The matrix elements of A and K are related 



Kjk = i [A, [A, A]].^ = i{xj - XkfAjk (5J 

17 



Finally, solving ( |5BD for Aj^ and putting into (0) we obtain 

4 l^j \ ' J/ I 

The first two terms are kinetic, coming from M^, while the last one is potential. 
Therefore, the hamiltonian H is 

H = j:lp^+'^T.7§^ + T.Vix.) (60) 

Note that the eigenvalues are kinematically coupled by an inverse-square type poten- 
tial with the angular momentum degrees of freedom. The connection of the matrix 
model to the Calogero model along the lines presented here and below was first es- 
tablished in [^]. Also, the hamiltonian ( |60D has been proposed independently of the 



matrix model as an S'f/(A^)-generalization of the classical Calogero system ||3g, |37 
We can now examine special cases: 

a) The most 'gauge invariant' sector is, of course, the one in which the angular 
momentum charges vanish, that is, J = 0. In that case, ( pOD for K = {] becomes the 
hamiltonian of non- interacting particles in an external potential V{x). This would 
be the case of 'standard' particles. 

b) For the next simplest case the angular momentum J should be as trivial 
as possible without vanishing. Only its eigenvalues are really relevant, since we 
can always perform a time-independent unitary transformation V which would shift 
U -^ VU and would rotate J -^ VJV~^ . The simplest choice would be to take the 
eigenvalues of J to be equal. Unfortunately, this is not possible since the traceless 
condition would make them vanish. The simplest possible choice is to take all the 
eigenvalues equal to £ except one, which would cancel the trace by being (1 — N)L 
This can be written in terms of an arbitrary column A^- vector v as 

J = i{vv^ - 1) , vh = N (61) 

in which case K becomes 

K = e{uu^ - 1) , u = U-\ (62) 

From (|58| ) we see that Ka = (no sum on i) and thus 

UiU* = 1 (no sum) (63) 

So the coefficient of the inverse-square potential in ( |60D becomes 

Finally, ( |60|) becomes 



KijKji = luiU* lujU* = f (z ^ j) (64) 
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This is the Calogero model! The potential strength g = P is related to the conserved 
charge i. Quantum mechanically, picking this charge will amount to a choice of 
statistics. We also get, at this stage, an arbitrary external potential V{x). 

c) More general choices of J amount to more variety in its eigenvalues. KijKji 
now, in general, becomes time- dependent and the dynamics more complicated. We 
postpone the discussion for the quantum case where it will be shown that this 
corresponds to Calogero particles having also internal degrees of freedom. This 
will be a generalization of the discussion of the first section, with irreps of SU{N) 
substituting the irreps of S]y. 

Now that we have this new approach we can use matrix technology to demon- 
strate the integrability of the Calogero model [^5|, ^. For V{x) = the matrix 



motion becomes free and M is conserved. The conjugation-invariant quantities 

J„ = trM" (66) 

are also conserved and in involution (the matrix elements of M are momenta and 
have vanishing Poisson brackets). From (|53) and (pSf) we have 



M,fc = U-'(5,kXj-il-5,k)^^^]u (67) 

\ Xj Xfc J 

= U-'(6,,x,-{l-6,,)^^^]u (68) 

When the above expression is inserted in the trace /„ = trM" clearly U drops 
and products of the form UiU* Ujul . . .u* will appear which reduce to powers of 
i. Therefore, the /„ reduce to expressions involving only Xi, Xi and the coupling 
constant £. These are the conserved integrals of the Calogero model. 

Starting from the matrix model the actual motion of the Calogero model can be 
obtained. The solution for M is 

M = B + Ct (69) 

for arbitrary matrices B, C. The conserved charge becomes 

J = i[M, M] = i[B, C] = ii{uu^ - 1) (70) 

By unitary transformations we can choose the phases of u such that Mj = 1; choices 
for B, C, then, that satisfy (^) are 

a 

Bjk = Sjk Qj , Cjk = SjkPj — (1 — 5jk) (71) 

Qj - Qk 

Qi and Pi are the initial conditions for Xi and Xi at time t = 0. Diagonalizing, then, 
M = B + Ct for the above B, C produces the motion of the system. Another choice 

is 

it 
Bjk = ^jkO-j + (1 — ^jk)-p, ^ ! Cjk = SjkPj (72) 
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Pi and ttj are asymptotic momenta and impact parameters. For t -^ ±00 the off- 
diagonal elements of B produce a perturbation of order t~^ to the eigenvalues, so 
the motion is determined by the diagonal elements Oj + Pit alone. We recover the 
result that the motion at asymptotic regions is the same as if the particles were free. 
We conclude by mentioning that the matrix model is also integrable and solvable 
in the presence of a harmonic oscillator potential V{x) = \uj'^x'^ . The non-hermitian 
matrix Q = M + iujM evolves as 

Q{t) = e^"*g(0) (73) 

and the matrix Q^Q is conserved. We leave it as an exercise to derive the conserved 
integrals and the motion of the corresponding Calogero problem. 

External potentials with up to quartic dependence on x also lead to integrable, 
although not so solvable, models [^. It is an open question to prove this is all there 
is, or to find yet more integrable potentials. For the purposes of statistics this is a 
rather academic issue. 

Finally, we may wonder what restricts us to one dimension. We chose a model 
with one matrix, and its eigenvalues corresponded to coordinates of particles on the 
line. We could, indeed, start with an appropriate model with many matrices, which 
would reproduce particle motion in higher dimensions [^|. The integrability and 



solvability properties of such extended models, however, are much less pleasant. The 
question of whether they represent a workable extension of statistics remains open. 



4.2 The unitary matrix model 

The hermitian matrix model works well for particles on the line but has trouble 
representing particles on periodic spaces. The most natural candidate for such 
models would be a unitary N xN matrix U. Its eigenvalues are phases and naturally 
live on the circle. We start, therefore, with a lagrangian that represents the invariant 
kinetic energy on the space of such matrices: 

C = -ItriU-'iir (74) 

A potential could in principle be included but we are interested in the translationally 
invariant case and will omit it. The treatment is similar as before, and we just 
summarize the relevant facts. 

The lagrangian is, in fact, invariant under separate left- and right-multiplications 
of U by time-independent unitary matrices and and there are two corresponding 
conserved matrix angular momenta L and R: 

U^VU: L = lUU-^ (75) 

U -^ UW-^ ■ R = -iU-^U (76) 

20 



The unitary conjugation that preserves the eigenvalues corresponds to W = V and 
its generator is 

J = L + R = i[U,U-^] (77) 

The rest of the discussion is as previously. Parametrizing 

U = VAV-'^ with A = rfmc/{e^^%...e^^^} (78) 

the hamiltonian becomes, after a few steps, 

^ = E^p^ + ^Er^fe- (79) 

where, as before, 

K = V-^JV (80) 

Choosing J = K = reproduces free particles on the circle, while choosing J = 
i{uu^ — 1) we obtain KijKji = (? and we recover the Sutherland inverse-sine-square 
model 

1 1 /2 

iJ = V-i,' + -V :—^ (81) 

This model is integrable and solvable by the same techniques as the hermitian 
one. The conserved invariant quantities are 

/„ = trL" = ii{-RY = iiiiU^^Uy (82) 

and the solution is 

U = Be'^' (83) 

with B a unitary and C a hermitian matrix satisfying 

BCB-^ -C = J (84) 

For the Sutherland case with J = (i{uu^ — 1), Ui = 1, B,C become 

a 

Bjk = Sjke'"' , Cjk = SjkPj + (1 - ^ifc) gifa-g,)_^ (85) 

where, clearly, qi and Pi are initial positions and momenta. 

We conclude by mentioning that, upon scaling x -^ ax, t -^ aH, the Sutherland 
model goes over to the free Calogero model. This is the 'infinite volume' limit. 
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4.3 Quantization and spectrum 

We will, now, perform the quantization of this system on a periodic space, using 
the unitary matrix model. We begin by defining a canonical momentum matrix 
conjugate to the 'coordinate' U 

n = ^ = -U-^UU-^ (86) 

dU 

The Poisson brackets are 

{f/j-fc, H/m} = ^jm^lk (87) 

n is somewhat unpleasant, being neither unitary nor hermitian. We prefer to work 
in terms of the hermitian matrices L and R defined previously 

L = iUU-^ = -iUU , R = -lU-^U = mu (88) 

Using (^) we derive the following Poisson brackets: 

{Ljk, Lim} = i{Ljm6ik — 6jmLik) (89) 

{Ljk^Rlm} = (90) 

{Rjk,Rlm} = iiRjmSlk — SjmRlk) (91) 

The above is nothing but two copies of the U{N) algebra in disguise. To see this, 
expand the matrices L and R in the basis of the fundamental generators of SU{N) 
T"- plus the unit matrix: 

L = L° + 2 J2 L^T" (92) 

a=l 

R = R° + 2 J2 R^T" (93) 

a=l 

with L°, L"', R", i?" numbers. Then use the SU{N) commutation relations 

as well as the normalization 

tr(T'^T^) = ^6ab (95) 

to show that the expansion coefficients satisfy the Poisson algebra 

{L^L^} = r^^L'^ (96) 

{L",i?^} = (97) 

{R^,R''} = r^^R" (98) 

while L°, R° are central. Note that the U{1) charges 

L° = -R° = ltr(-zf/-if7) = ^ E i^ (99) 

i 
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are essentially the total momentum of the system. 

We are now ready to perform quantization. In the [/-representation, where states 
are functions of U, U becomes the matrix derivative Ujk = —iSu, acting as 

6utT{UB) = B , 6utT{U-^B) = -U'^BU'^ (100) 

where i? is a constant matrix, and similarly on expressions containing more f/'s. L 
and R, upon proper ordering, are represented as 

L = -U6u , R = 6u ■ U (101) 

where in R it is understood that we first act with the derivative and then right- 
multiply the result by U. With this ordering, L and R become the proper U{N) 
operators acting as 

Ltr{UB) = -UB , Ltr{U-^B) = BU~^ (102) 

RiiiUB) = BU , Rti{U-^B) = -U'^B (103) 

It is also useful to express their action on arbitrary functions of U as 

tr(*eL)/(f/) =/((l-^e)f/)-/(f/) (104) 

tT{ieR)f{U) = f{U{l + te)) - /(t/) (105) 

where e is an arbitrary infinitesimal hermitian matrix, emphasizing their role as 
generators of left- and right-multiplication on U. Correspondingly, the operators L" 
and i?" satisfy the SU{N) algebra. Their action can be obtained by taking e = eT" 
with e an infinitesimal scalar parameter, that is, 

tsL'^fiU) =f[{l-reT-)U)-f{U) (106) 

leR'^fiU) =f{U{l + ieT-))-f{U) (107) 

The hamiltonian, being classically the kinetic term on the manifold of unitary ma- 
trices U{N), quantum mechanically becomes the laplacian operator on the manifold 
50]. Using (P^IPBD it is expressed as 



H = itrL^ = Y^iLr + lN{Lr = ^(i?^)^ + ^{R^ = hrR' (108) 

It is, therefore, the common quadratic Casimir of the left- and ngh.t-SU{N) algebra 
plus the square of the U{1) charge, the two parts identifiable as the relative and 
center-of-mass energy respectively. 

Quantum mechanical states grouping into irreducible representations of the L 
and R SU{N) algebras will, thus, be degenerate multiplets of the hamiltonian. The 
U{1) (center of mass) part trivially separates: we can boost any state by any desired 
total momentum NP by multiplying the wavefunction by (det U)^ . We will examine 
only the SU{N) part from now on. 
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A natural basis of states for the Hilbert space are the matrix elements of the 
unitary irreducible representations (irreps) of SU{N). Let R denote such an irrep, 
R{U) the matrix that represents U in this irrep and RapiU) the a(5 matrix element 
of this matrix. Clearly a and (3 range from 1 to the dimensionality of i?, d^. Rap{U) 
are a complete orthonormal basis of wavefunctions for [/, that is 

[dU]R^p{U)R'^,{Uy = 5rr,5^^5^s (109) 

with [dU] the volume element on the space of SU{N) matrices as implied by the 
metric ds'^ = —ti{U~^dU)'^, also called the Haar measure. 

We will, now, show that each R^piU) is an eigenstate of the hamiltonian with 
eigenvalue equal to the quadratic Casimir of -R, Cr. Qualitatively, after the discus- 
sion of the last paragraphs, this should be obvious: L and R generate the trans- 
formations U -^ V'^U and U -^ UW. R{U) transforms in the conjugate irrep R 
under L and in the irrep R under R. Since H is the common quadratic Casimir of L 
and R we conclude that all dj^ states RapiU) are energy eigenstates with eigenvalue 

Cr = Cr. 

(If you are confused about L generating U -^ V~^U rather than U —>■ VU, think 
of the difference between active and passive transformations, which is relevant when 
shifting from classical to quantum: iIj{x — a) shifts the wavefunction by +a. Also, 
although classical transformations on U compose properly, 

Vi{V2U) = {V^V2)U (110) 

quantum mechanically the operators V that perform the shift U -^ VU on the 
argument of the wavefunction would compose 

Vi{V2f{U)) = VifiV^U) = fiV^V^U) = {V2V1) f{U) (111) 

Therefore we need to invert the action of V to get the right composition law. See 
also the discussion of irreps of Sn in the section about permutation group statistics.) 
Let us prove the fact HRap{U) = CRRapiJJ) more analytically. Since R{U) is a 
representation, it obeys the group property 



Rap{UV)=Y.R^^R^p{y) (112) 



From ( |106| ) we have 



(1 + ie)L''R^p{U) = Ra/siil - ieT'')U) = Ra^{l - ieT'')R^p{U) (113) 
= R^p{U) - ieRl^R^p{U) (114) 

where i?" = R{T°') is the a-th generator of SU{N) in the R representation. So 

L'^R^piU) = -Rl^R,p{U) (115) 
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and 

^(L«)2i?,^(f/) =J2R^^^R;,RspiU) = J2iR^)lsRsfsiU) (116) 

a a a 

The sum I]a(-R")^ appearing above is the quadratic Casimir in the irrep R and is 
proportional to the identity matrix 6a5- So, finally, 

HR^^iU) = CrR^p{U) (117) 

Incidentally, the spectrum spanned by Cpt for all R is nothing but the spectrum of 
N free fermions on the circle with the ground state energy and the center-of-mass 
energy subtracted, where the lengths Ri of the rows of the Young tableau of R 
correspond to the "bosonized" fermion momenta 

Pi = Ri-t + l (118) 

and where the center-of-mass energy has been subtracted. The condition Ri > 
Ri+i for the rows amounts to the fermionic condition pi > Pi+i. The spectrum of 
the full matrix model, then, is identical to the free fermion one but with different 
degeneracies. 

We have, therefore, identified all energy eigenstates of the matrix model. It 
remains to implement the quantum analog of the choice of angular momentum J, 
identify the corresponding reduced quantum model, and pick the subspace of states 
of the full model that belongs to the reduced model. 

J obeys itself the SU{N) algebra (it is traceless, no U{1) charge). A choice of 
value for J amounts to a choice of irrep r for this algebra. States within the same 
irrep are related by unitary transformations of U and give the same dynamics; they 
will be discarded as gauge copies, and only the choice of irrep will be relevant. Since 
J = L + R, we see that states transforming under (L, R) in the {R, R) irreps will 
transform in the R x R under J. So, only irreps r that are contained in the direct 
product of two mutually conjugate irreps can be obtained for J. This amounts to 
irreps r with a number of boxes in their Young tableau that is an integer multiple 
of N. (To get a feeling of this, consider the case N = 2. Then J is an orbital-like 
realization of the angular momentum through derivatives of U and clearly cannot 
admit spinor representations.) 

We must, therefore, project the dj^ states in RajsiU) to the subspace of states 
transforming as r under L+R. Call G{R, a; R, f3\r, 7) the Clebsch-Gordan coefficient 
that projects these states to the 7 state of r. Then the relevant states for this model 
become 

^RiU) = J2 Raf3{U)G{R, a; R, (3\r, 7) (119) 

The index 7 labeling the states within r, as we argued before, counts the dr gauge 
copies and does not imply a true degeneracy of states. The degeneracy of the states 
produced by each R is, then, given by the number of times that the irrep r is 
contained in the direct product R x R or, equivalently, the number of times that R 
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is contained in R x r. Calling this integer D[R,r] R), we obtain for the spectrum 
and degeneracies: 

Er = Cr, DR = D{R,r;R) (120) 

In particular, if Dr = the corresponding energy level is absent from the spectrum. 

Concluding, we mention that an approach which also reproduces the spectrum 

and states of the Sutherland model is two-dimensional Yang-Mills theory on the 

circle PTl O. This approach is essentially equivalent to the matrix model above 



and we will not be concerned with it. 

4.4 Reduction to spin-particle systems 

So we have derived the spectrum, degeneracy and wavefunctions of the matrix model 
restricted to the sector J = r. Classically these restrictions represented free particles 
(J = 0), Sutherland particles (J = £(ft'^ — 1)) or something more general. What 
are the corresponding quantum systems? 

To find these, let us reproduce here the expression of the reduced hamiltonian 
in one of these sectors: 

^-Eirf + iEjIpfe-s. (121) 

i ^ ^ i^j ^ sm ^2~ 

This expression remains valid quantum mechanically upon a proper definition (or- 
dering) of the operator K. The only residual quantum effect is a constant term 
Eo that comes from the change of measure from the matrix space to the space of 
eigenvalues. 

Let us expand a bit on this without entering too deeply into the calculations. 



(For details se, e.g., @].) The Haar measure in terms of the diagonal and angular 
part of f/ has the form 

[dV] = A^ [dV] (122) 

where [dV] is the Haar measure of V and A is the Vandermonde determinant 

A = n2sin^^^^ (123) 

i<j 

To see this, write the 'line element' — tr(t/^^c?t/)^ in terms of V and Xi using (|78D 
and obtain 

- tT{U-^dUy = ^ t/x^ - ^ 4 sin^ ^^i^^(V''dV)ijiV-'dV)ji (124) 

This metric is diagonal in dxi and (y~^dV)ij. The square root of the determinant 
of this metric, which gives the measure (volume element) on the space, is clearly A^ 
times the part coming from V which is the standard Haar measure for V. (We get 
two powers of 4sin^ x^-xj ^^ ^^^ determinant, one from the real and one from the 
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imaginary part of {V ^dV)ij, so the square root of the determinant has one power 
of A2.) 

To bring the kinetic x^-part into a 'flat' form (plain second derivatives in Xi) 
we must multiply the wavefunction with the square root of the relevant measure 
(compare with the change from cartesian to radial coordinate in spherical potential 
problems). The net result is that the wavefunction \E' in terms of Xi and V is the 
original wavefunction ip{U) of the matrix model times the Vandermonde determi- 
nant. This, however, also produces an additive constant Eg which comes from the 
action of the entire Xj-kinetic operator on A. Noticing that A is nothing but the 
ground state wavefunction of A^ free fermions on the circle, we see that Eo is the 
relevant fermionic ground state energy 

Eo = ^ ^^ ' (125) 

This is the famous 'fermionization' of the eigenvalues produced by the matrix model 
measure. 

To determine the proper ordering for K we examine its properties as a gener- 
ator of transformations. Since U = VAV~^, and J generates U -^ V'UV'~^ = 
iy'V)A(y'V)~^, we see that J generates left-multiplications of the angular part V 
oi U. K = V~^JV, on the other hand, generates n5f/i/;-multiplications of V, as can 
be see from its form or by explicit calculation through its Poisson brackets. As a 
result, it also obeys the SU{N) algebra. Its proper quantum definition, then, is such 
that it satisfies, as an operator, the SU{N) algebra. It clearly commutes with the 
diagonal part Xi and its momentum pi, since it has no action on it. Its dynamics are 
fully determined by the hamiltonian ( |121D and its SU{N) commutation relations. 



We can, therefore, in the context of the particle model (|121| ), forget where K 
came from and consider it as an independent set of dynamical SU{N) operators. K, 
however, obeys some constraints. The first is that, as is obvious from K = V~^JV, 
K carries the same irrep r as J. The second is subtler: a right-multiplication of 
V with a diagonal matrix will clearly leave U = VAV~'^ invariant. Therefore, this 
change of V has no counterpart on the 'physical' degrees of freedom of the model 
and is a gauge transformation. As a result, we get the 'Gauss' law' that physical 
states should remain invariant under such transformations. Since K generates right- 
multiplications of V, and Ku (no sum) generates the diagonal ones, we finally obtain 

(no sum) Kii = (on physical states) (126) 

(A more pedestrian but less illuminating way to see it is: J = i[U^^,U], being a 
commutator, vanishes when sandwiched between the same eigenstate of U. Since K 
is essentially J in the basis where U is diagonal, its diagonal elements vanish.) Note 
that the constraint ( |126|) is preserved by the hamiltonian (|121|) . 



The above fully fixes the reduced model Hilbert space as the product of the 
A'^-particle Hilbert space times the (i,.-dimensional space of K, with the constraint 
( |126| ) also imposed. The further casting of the model into something with a more 
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direct physical interpretation relies upon a convenient realization of K. Any such 
realization will do: simply break the representation of SU{N) that it carries into 
irreps r and read off the spectrum for each r from the results of the previous section. 
We shall implement K m. & construction a la Jordan-Wigner. Let a^i, al^i, 
771 = 1, ... g, -i = 1, ... A^ be a set of Nq independent bosonic oscillators [Q: 

Then 

K''=J2 (^miTt.a^j (128) 

m=l 

is a realization of the SU{N) algebra. (T^°- are the matrix elements of T"-.) The 
corresponding matrix elements of K are 

^ij = Yl \ (^Irnamj - ^ ( I] aL«mfc j 5ij \ (129) 



Correspondingly, the coefficient of the Sutherland potential in (|121|) is (for i ^ j) 



m,n 



We already see that the degrees of freedom of K are redistributed into degrees of 
freedom for each particle in the above. Specifically, defining 

1 f '^ \ 



q 



\s=l 



and comparing with (|129|) we see that the Si are A^ independent sets of operators 
each satisfying the SU{q) algebra. Before expressing KijKjt in terms of the Si let 
us see what the constraint (|126|) implies: 

m=l m,k 

J2m,k (Knk'^mk commutcs with all Kij and all Si^mn- It is, therefore, a Casimir and can 

be chosen as a fixed integer iN equal to the total number operator of the subspace 

of the oscillator Fock space in which the model lives. The above constraint, then, 

implies 

q 

J2 alaiami = i (133) 

m=l 

(We see why we had to choose the total number operator to be a multiple of A^: 
the operator in ( |133| ) above is also a number operator and can have only integer 
eigenvalues.) Using this in ( |131| ) we can express 
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and therefore 

where Si ■ Si = ti^SiSj) is the S'f/(g)-invariant scalar product of the two SU{q) 
'vectors.' We finally obtain the hamiltonian as [E2[ 



g=E;p?+;E ,J..-.; (136) 

So it is a Sutherland-like model but where the particles also carry SU{q) internal 
degrees of freedom ('spins') and the potential contains a pairwise antiferromagnetic 
interaction between the spins. 

It remains to specify the representation in which the SU{q) spins are and find 
the irreps contained in this realization of K, therefore obtaining the spectrum. A 
realization of the form ( |129| ) for g = 1 in terms of bosonic oscillators contains all 
totally symmetric irreps of S{N) (that is, the ones with a single row in their Young 
tableau). (|129[) is essentially the direct product of q such independent realizations, 



so it contains all direct products of q totally symmetric irreps. This produces all 
irreps with up to q rows in their Young tableau, perhaps more than once each. 



The constraint ( |133|) , however, implies that the total number of boxes in the Young 
tableau of these irreps is iN. We recover once more the constraint that we derived 
before based on the origin of r as a component oi R x R. 

Similarly, the realization ( |131| ) of Si contains all the totally symmetric irreps of 
SU{q). ( |133|) implies that the number of boxes of these irreps is equal to i, so the 
spins Si are each in the £-fold symmetric irrep of SU{q). Solving this model amounts 
to decomposing the tensor product of these A^ spins into irreducible components of 
SU{q). Each such component corresponds to a subspace of the Hilbert space with 
a fixed total spin S. This same irrep, interpreted as an irrep r of SU{N), will be 
the corresponding irrep of K, and also of J, and thus will determine the spectrum 
of this sector through (|120| ) . 

Let us elucidate the above by reproducing the two simplest cases: free particles 
and (spinless) Sutherland particles, comparing with the classical treatment. 

a) Free particles correspond to J = K = 0. So there is no spin and no potential 



and we have non-interacting particles. From ( [12CI| ) we see that all Dr are one, and 
thus the spectrum is the free fermion one, as commented before. The matrix model 
naturally quantizes free particles as fermions. 

b) Spinless Sutherland particles correspond, classically, to J = i{vv^ ~ !)• So 
J is rank one (ignoring the trace). Quantum mechanically this corresponds to the 
irrep r of J having only one row and therefore only one independent Casimir. Since 
q in the realization above corresponds to the number of rows, we must have q = 1- 
Spins, therefore, are absent. The strength of the potential becomes i{i + 1) where 
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IN is the number of boxes in the one row of r. By standard Young tableau rules we 
see that the degeneracy Dr is one if the row lengths of R satisfy 

Ri>Ri+i + ^ (137) 

else it is zero. The spectrum of this model is, then, the same as the spectrum of free 
particles but with the selection rule for their momenta 

Pi>P^+l + l+l (138) 

We recover the 'minimum distance' selection rule of the CSM model that led to 
the interpretation as particles with generalized statistics! Only, in this case, the 
statistics parameter £ + 1 is a positive integer. 

We mention here that a Jordan- Wigner realization of K in terms oi fermionic os- 
cillators is also useful and leads to particles with spins interacting via ferromagnetic 



Sutherland- type potentials. The hamiltonian becomes [42 



H = y^^-p^-W ' ;: / (139) 

i -^ ^ i^j ^ sm -^— 

where now the spins are in the £-fold anteymmetric irrep of SU{q). We will not 
elaborate further and leave the details as an exercise to the reader. 

In conclusion, the matrix model has provided us with the following: 

1. An augmentation of the permutation group into the SU{N) group and a 
corresponding possibility to define statistics through the irreps of SU{N). 

2. A realization of generalized scalar statistics but with quantized statistics 
parameter £ + 1 in terms of the CSM model. 

3. A realization of generalized 'non-abelian statistics' in terms of particles with 
internal degrees of freedom interacting through a generalized CSM-type potential. 

4. A systematic way of solving the above models. 
What the matrix models has not provided is 

1. A realization of generalized statistics for fractional statistics parameter. 

2. A realization of spin-CSM systems with the spins in arbitrary (non-symmetric) 
representations. 

3. A control of the coupling strength of the potential for the spin-CSM models. 
(Note that the coefficient of Si ■ Sj terms is fixed to ±2 and also the constant therm 
i{i + q) is entirely fixed by the spin representation.) 

There exist generalizations of the above models both at the classical 



and the quantum level |^. They all share, however, the limitations spelled out 
above, especially (3). These restrictions are important, in the quest of more general 
statistics but also from the more practical point of view of solving spin-chain models 
with spins not in the fundamental representation, as we will shortly explain. For 
this reason, yet a different approach will be pursued in the next section, namely, the 
operator approach. 
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5 Operator approaches 

The matrix model connection provided us with a powerful tool that not only allowed 
us to generalize statistics but also led to the full quantum solution of a set of spin- 
generalized CSM models. 

As noted, however, in the conclusion of the preceding lecture, the matrix model 
fixes the coefficient of the spin-interaction and scalar interaction terms to ±2 and 
±i{i ± q) respectively. We cannot choose these coefficients at will. 

We would like to have an approach that defeats this restriction and leads to spin 
models with arbitrary coupling strengths. (This is necessary to attack spin-chain 
systems through the infinite-coupling limit trick to be explained later.) Such an 
approach should also be able to bypass the excursion to matrix models and deal 
more directly with these systems in an algebraic way. This will be achieved with 
the exchange operator formalism [^ (also known as the Dunkl operator formalism 



5.1 Exchange operator formalism 

Consider the operators Mij that permute the coordinate degrees of freedom of A^ 
particles in one dimension which could, in principle, also have internal degrees of 
freedom (M for metathesis, to avoid confusion with momenta Pi). They satisfy the 
permutation algebra (symmetric group), in particular 

M,, = Mr/ = Ml = M,, (140) 

[Mij, Mki] =0 if 2, j, k, I distinct (141) 

MijMjk = MikMij if i,i, k distinct (142) 

Any operator Ai on the phase space satisfying 

MijAk = AkMij if i,j,fc distinct (143) 

Mi.Ai = AjMij (144) 

will be called a one-particle operator (even though it may involve the coordinates 
and momenta of many particles). 



We construct the following one-particle operators [^ 



^ri=P^ + J2 ^^(^i - ^j)MiJ -P^ + T. i^ijM^j (145) 

We shall view the ttj as generalized momenta. To ensure their hermiticity the pre- 
potential W{x) should satisfy 

W{-x) = -W{xy (146) 
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We shall construct the corresponding 'free' hamiltonian from tTj 



H = T.U 



(147) 



In terms of the original pi this hamiltonian will, in general contain linear terms. To 
ensure that such terms are absent we must further impose 



W{-x) = -W{x) = real 



148) 



With the above restriction the hamiltonian H and commutation relations of the ttj 
become 

[vTi, TTj] = Y^ Wijk{Mijk - Mjik) 



1 



H=J:oPl+J:MJ+KM^,)+ E W^,'^M, 



i ^ i<j i<j<k 

where we defined the three-body potential and cyclic permutation 



Wijk 






(149) 
(150) 



(151) 
(152) 



To obtain an interesting and tractable model, Wijk-, which appears in the commuta- 
tor [ttj, Tij\ and also as a three-body potential, should vanish or at most be a constant. 
This leads to a functional equation for W{x): 



W{x)W{y) - W{x + y) [W{x) + W{y)] = const (= Wijk) 



(153) 



We present the solutions: 

a) Wijk = -^ W{x) = i/x 

b) Wijk = -f <0 -^ wlx) 

c) Wijk = +f>0 -> W{x) 
Let's examine each case. 



?cot a; 
' coth X 



a) In this case the ttj become [E3, H^, 



TTi 



Pi + Y. —M, 



j^i ^ij 



U 



and satisfy 



[tTj, Tij] = 



(154) 



(155) 



The TTj commute, so we can consider them as independent momenta. The hamilto- 
nian reads 

1 , ^i{i-Mi/ 



H = Y^P^ + Y 



^ijj 



Kj 



Xf 



(156) 
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We obtain a Calogero-like model with exchange interactions. Yet it is nothing 
but a free model in the commuting momenta ttj. Integrability is immediate: the 
permutation-invariant quantities 



E 



vr,- 



(157) 



obviously commute with each other. If we assume that the particles carry no internal 
degrees of freedom and are bosons or fermions then Mij = ±1 on physical states. The 
model becomes the standard Calogero model and we have proved its integrability 
in one scoop. (You may be left with a question mark: the hamiltonian and the 
other integrals J„ become the standard Calogero ones if M^j = ±1, so these reduced 
integrals will commute on the bosonic or fermionic subspace; but will they also 
commute on the full Hilbert space? Prove for yourself that this is indeed the case.) 



We can also construct harmonic oscillator operators [47, 4£]. The commutators 



between Xi and ttj are 



Xi, 71 i -- 


= I l+^Y.M^J 




\ i^» 


Xi, 7Tj -- 


= -im,, {1^3 



(158) 
(159) 



Defining 



72 
1 

71 



[11 j — lUX,: 



(TTi + iuJXi 



(160) 
(161) 



we can show 



[a„al] = u(l + iY.M^j 

[ai,a]] = -uiMij {ij^j 
[ai,aj] = [a|,a]] = 



(162) 

(163) 
(164) 



This is an extended version of the Heisenberg algebra involving the permutation 
operators. The corresponding oscillator hamiltonian reads 



H = E 



1 



1 



1 



a tti 



+ Oial) = ^ -p1 + ^ -uj'^Xi + ^ 



M,,) 



i<j ^ij 



(165) 



and satisfies 



[if, ai] = uai , [if, aj] = uja\ 



(166) 
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This is the harmonic Calogero model with exchange interactions, which becomes 
again the standard model on bosonic or fermionic subspaces for particles without 
internal degrees of freedom. Since 

H = J2 (^U + \nuj + Uu Y. M,, (167) 

we see that on bosonic or fermionic spaces the state annihilated by all Oj (if it 
exists) will be the ground state. Solving a^^ = we obtain for the ground state 
wavefunction 



i^B = ni^.,iv^-s^^>' (168) 

i<j 

i^F = Y[{^gn{x,,)\x,,\-'}e-"^'^^^^" (169) 

i<j 

For i > the bosonic state is acceptable, while for £ < the fermionic one is 
acceptable. In the "wrong" combinations of statistics and sign of i the ground state 
is not annihilated by the Oj, but it is still annihilated by all permutation-invariant 
combinations of the aj. 

From (|166|) we see that we can find the spectrum of this model for fermions 
or bosons by acting on the ground state with all possible permutation-symmetric 
homogeneous polynomials in the a]. A basis bor these is, e.g.. 

An = E(4)" (170) 

i 

So the spectrum is identical to non-interacting fermions or bosons, but with a dif- 
ferent ground state energy. For the 'right' combinations of i and statistics, where 
( |169| ) are the correct ground state wavefunctions, the ground state energy is 

which is the correct Calogero result. 
Finally, the quantities 

4 = E^r = E(4«.r (172) 



can be shown to commute [47|, and therefore this system is also integrable. It is left 
as an exercise to find the commutation relations of the hi and show that [In, Im] = 0. 
b) In the case W{x) = icotx we have 

TTj = pi + i cot XijMij (173) 

K, vr,] = -f E(^ufc - M,.k) (174) 

k 
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so the momenta are now coupled. The hamiltonian becomes 

. z, .^j bin Xij y z i<j<k I 

We obtain the Sutherland model with exchange interactions plus an extra term. 
On bosonic or fermionic states this becomes an overall constant and we recover the 
standard Sutherland model. Again, since H is by construction positive definite, if a 
state satisfying ttj?/' = exists it will be the ground state. We obtain 

^F = Ylsgn{xij)\ sin Xij f (177) 

i<j 

which are acceptable for the same combinations of i and statistics as before. For 
both cases Mijk = 1 so 

E. ^ f^^^^ (178, 

is the correct Sutherland model ground state energy. The excited states can again 
be obtained in a (rather complicated) algebraic way |5D|. Finally, the quantities 



2/ 
vf, = TT, + ^ E M,, =p, + e-- E ;:i^r— -M,, (179) 

can be shown to have the same commutation relations as the hi defined previously 
for the harmonic system. Therefore, the integrals constructed from them 

/n = E^r (180) 

i 

commute and the model is integrable. 

c) For W{x) = £ cothx we have a similar commutation relation and a hamiltonian 

h^e^.j: ^(i^ . f f ^^^ . E mJ ,181) 

i 2 i^j smh Xi, \ 2 .^^.^^ J 

This is the inverse-hyperbolic-sine-square model and supports only scattering states. 
Its integrability can be obtained as for the Sutherland model above, or simply as 
an 'analytic continuation' of that model for imaginary period of space. We will not 
examine it any further. 

In conclusion, an exchange-family of models was introduced, solved and related 
to the standard CSM models in spaces of definite symmetry. It is remarkable that 
all these proofs work directly, and only, at the quantum domain (there is no classical 
analog of Mij). 
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5.2 Systems with internal degrees of freedom 

We can easily extend the previous results for particles with internal degrees of free- 
dom. For this, assume that the particles are distinguishable or, equivalently, that 
they carry a number q of (discrete) internal degrees of freedom (species) that can 
be used to (partially) distinguish them. Their states are spanned by \x,a >, where 
cr = 1, . . . g counts internal states. The total permutation operator Tij, then is 

Tij = Mijaij (182) 

where aij is the operator that permutes the internal states of particles i and j. 

Let us, then, simply take states that are bosonic or fermionic under total particle 
exchange: Tij = ±1. On such states 

Mij = ±(Tij (183) 



and the Calogero and Sutherland exchange model hamiltonians become pT 



We get the Calogero and Sutherland models with spin-exchange interactions. From 
the completeness relation for the fundamental SU{q) generators T" 

9'-i ^ „ 1 1 

E '^apT^S = -^^asSyp - —Sap6^5 (186) 

0=1 ^ ^y 

we deduce the form of the operators cTij 

aij = 2Si ■Sj + - (187) 

where 5"" acts as T" on the internal states of particle i. So the spin-dependent 
interaction coefficient of the potential in the hamiltonian takes the form ^, ^ |5l|. 
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^il2Si-SjTi+-] (188) 

We have recovered the ferromagnetic and antiferromagnetic spin model of the pre- 
vious section but with arbitrary coefficient! On the other hand, the spins are nec- 
essarily in the fundamental of SU{q). So we have obtained a generalization of the 
coupling constant with respect to the matrix model but a restriction of the allowed 
spins. 

36 



Note that d, here is an arbitrary parameter, while £ in (|136|) was the size of the 



symmetric representation of Si. For i = 1 and spins in the fundamental, the matrix 
model and exchange-operator model agree. It is interesting to note that we can go 
from ferromagnetic to antiferromagnetic interactions either by changing the sign of 
a or by changing the statistics of the particles. 

The solution of the above models can be obtained algebraically. For the spin- 
Sutherland model this is rather complicated and is related to the so-called Yangian 



symmetry |55, pq, 5/, 58|. For the spin-Calogero model it is easier [M]. Let us 



concentrate on the model with interaction £(— 2S'j-S'j=f£— -) and define the operators 

i i 

and their hermitian conjugates. These form a complete set for all permutation- 
symmetric creation and annihilation operators for all species of particles. Yet the 
commutators among themselves and with H do not involve i. They create, therefore, 
the same spectrum of excitations over the ground state as A^ noninteracting bosons 
or fermions with q species. For i > the ground state is the bosonic one: 

i^B = n \x^/e-"^'^^^''hsiW^}) (190) 

i<j 

where Xs is a totally symmetric state in the ai. The set of all Xs forms the A^- 
fold symmetric irrep of the total spin S = J^iSi- Therefore the ground state is 
(A^ + q — 1)\/N\{q — 1)! times degenerate. For i < the above is not normalizable 
any more. But we remember that we can obtain the same model by starting from 
fermions and the opposite coupling —i>0. The ground state, then, is of a fermionic 
type 

^F = U-^f (m^.,.) n \x.,\-'^-^^e^-^^^^ (191) 

P \ i / i<j 

where P are total particle permutations and «« are a set of fixed values for the 
indices ctj that determine the state. Clearly the ground state will be obtained for 
the minimal total power of Xi appearing above, and that will happen for a maximally 
different set of values ««. These states form the n-fold antisymmetric irrep of the 
total spin S, where n = N{modq). The ground state is, thus, q\/n\{q — n)\ times 
degenerate. The above spectra will come handy later. 

5.3 Asymptotic Bethe Ansatz approach 

We already mentioned that there are elaborate algebraic approaches to derive the 
spectrum of the spin-Sutherland model, based on the Yangian symmetry. We will, 
instead, take a lower-key approach which reproduces the same spectra and is phys- 
ically more lucid, although not as rigorous. We will take the ABA route. 

Consider distinguishable particles of the exchange-Calogero type without external 
potential, coming in with asymptotic momenta ki and scattering off each other. 
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Before scattering, their positions are in some definite ordering determined by the 
ordering of their momenta (it is the inverse of that ordering). 

The key observation is that, after scattering, the particles have simply 'gone 
through' each other with no backscattering [Q. The impenetrable 1/x^ potential 
has become completely penetrable in the presence of the exchange term! You can 
prove this fact by examining the asymptotic properties of a simultaneous eigenstate 
of TTi, . . . TTjv which is obviously an energy eigenstate: at Xj -^ ±00 the prepotential 
terms are vanishing and we simply have eigenstates of the individual pi. Since 
there are no pieces with the values of pi permuted (coming from backscattering) 
we have complete transmission. (To explicitly see how it works, it is instructive 
to consider the two-body problem, decompose it into symmetric and antisymmetric 
parts, scatter and recombine the parts after scattering. A relative phase of vr between 
the two parts is what produces the effect.) 

(Puzzle: what happens with the correspondence principle? With h back in, the 
interaction coefficient is £(£ — TiMij). How can a term of order h produce such a 
dramatic effect, particles going through each other, in the /i — ^ limit?) 

So the only effect of the scattering is a phase shift of the wavefunction which, as 
we have said, is the sum of two-body phases 

This is true on an infinite space. On a periodic space we can still use the above 
result, together with the requirement for periodicity for the wavefunction, to derive 
the spectrum. This is the ABA method and is expected to reproduce the correct 
results in the thermodynamic limit of many particles at constant density 16^]. It 
gives, in fact, the exact answer for the Sutherland model ||26|, so we can expect it to 



work also in the present case. For a space of period 27r the result is 

2'Kki + ^^£sgn{ki — kj) = 2imi (193) 

■K 

The left hand side counts the total phase picked up by a particle going round the 

space and scattering off to the other particles in the way. Ui are arbitrary integers, 

ensuring periodicity. There are, however, some constraints on the choice of of ni 

that are imposed by continuity from the £ = case: 

-If ki < kj then Ui < Uj 

-If Ui = Uj there is a unique solution, that is, ki < kj and ki > kj represent the 

same state. 

These rules are important to avoid overcounting and to discard spurious solutions. 

With these, the spectrum obtained is the same as the one derived with more rigorous 

methods. For the ordering ni < . . . njsi the solution for ki is 

N + 1 
h = ni + l{i —) (194) 

38 



and similarly for other orderings. We see that the ABA momenta /cj are the same 
as the quasimomenta that we have previously defined. 

The bottom line is that the spectrum and degeneracies are the same as those 
of distinguishable particles obeying generalized selection rules for their momentum. 
Still, what fixes the degeneracy of states is the different ways that we can distribute 
the particles to the quantum numbers rii, rather than ki (see, especially, the second 
rule above). A state of A^ particles with the same rii, for instance, is nondegener- 
ate although they, seemingly, have different k^ which would imply a permutation 
degeneracy. 

For particles with spin the construction above, in combination with the trick of 
the previous subsection of starting with fermions or bosons, produces a spectrum 
with degeneracies the same as those of free particles (the rii are 'free' quantum 
numbers). As argued before, for ferromagnetic interactions we must choose bosons 
and combine their spins accordingly, while for antiferromagnetic interactions we 
must choose fermions. To spell it out, this means the following: 

1. Choose a set of quantum numbers rij. The ordering is immaterial, since we have 
identical particles, so you can choose rii < . . . tin- 

2. Place your particles on these quantum numbers and put their spins in the ap- 
propriate state. For the ferromagnetic case treat them as bosons: the total spin 
of particles with the same n, transforms in the symmetric tensor product of their 
spins. For the antiferromagnetic case treat them as fermions: the total spin of of 
particles with the same rii transforms in the antisymmetric tensor product of their 
spins; clearly up to q can have the same rii in this case. 

3. Calculate the energy of this state in terms of the ABA momenta ( |194| ): E = J2i kf- 
It should be obvious that similar rules applied to the spin-Calogero system reproduce 
the spectrum derived in the last subsection. This method can be used to calculate 
both the statistical mechanics (large N) of these systems and the few-body spectra. 

5.4 The freezing trick and spin models 

Now that we have a tractable way of solving spin-CSM systems with arbitrary 



strength of interaction we can introduce the freezing trick [35| and deal with spin 
chain models. 

Consider, first, the previous ferromagnetic or antiferromagnetic spin-Sutherland 
model. Take the limit i ^ oo. The potential between the particles goes to infinity, 
so for any finite-energy state the particles will be nearly 'frozen' to their classical 
equilibrium positions. In fact, even the excitation energies around that configuration 
will go to infinity: the ground state energy scales like P, while the excitations scale 
like Nin + n^ with n some excitation parameter. So, to leading order in i the 
spectrum becomes linear and of order i. These excitations correspond, essentially, to 
phonon modes of small oscillations around the equilibrium positions of particles. The 
'stiffness' of oscillations is, of course, proportional to the strength of the potential 
£^ and the spectrum is proportional to the frequency, of order i. 
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The quantum fluctuations of tlie particle positions in any state will scale like 
the inverse square root of the oscillator frequency, that is, like 1/Vt. But, in the 
hamiltonian, the piece coupling the spins to the kinematical degrees of freedom is 
proportional to 1/sin^Xjj. In the large-£ limit, thus, this term becomes a constant 
equal to its classical equilibrium value; so, in that limit, spin and kinematical degrees 
of freedom decouple. (Note that the spin part is also of order L) The hamiltonian 
becomes 

H = Hs + iH,p,n (195) 

with Hs the spinless Sutherland hamiltonian and Hspm the spin part 

H.,.n = TY.J^^ (196) 

where the classical equilibrium positions Xj are equidistant points on the circle: 

^3 = ^ (197) 

The hamiltonian ( |196| ) above describes a spin chain consisting of a regular peri- 
odic lattice of spins in the fundamental of SU{q) coupled through mutual ferro- or 
antiferromagnetic interactions of strength inversely proportional to their chord dis- 
tance. It is the well known SU{q) Haldane-Shastry (HS) model |61|, |62[. According 



to the above, its spectrum can be found by taking the full spectrum of the corre- 
sponding spin-Sutherland model in the large-£ limit, 'modding out' the spectrum of 
the spinless model and rescaling by a factor 1/i. Each state will inherit the spin 
representation of its 'parent' spin-Sutherland state. So, both the energy and the 
total spin of the states of the HS model can be determined this way. Commuting 
integrals of this model [Q can also be obtained this way [B3|. At the level of the 
partition function at some temperature T we have 

Z(£T) 
ZspUT) = lim -^ (198) 



From this, the thermodynamics of the spin chain model can be extracted ^0| . 

We will not give the details of this construction here. We urge anyone interested 
to solve this way a few-site (two or three) spin chain, see how it works and deduce 
the 'construction rules' for the spectrum of a general spin chain. Let us simply 
state that the many degeneracies of the spectrum of the HS model (larger that the 
total spin SU{q) symmetry would imply), which is algebraically explained by the 
existence of the Yangian symmetry, can, in this approach, be explained in terms of 
the degeneracies of free particles. (The degeneracies are not identical, due to the 
modding procedure, but related.) 

For the spin-Calogero model a similar limit can be taken, scaling also the external 
oscillator frequency as u; — > £a; to keep the system bound. The classical equilibrium 
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positions of this model are at the roots of the A^-th Hermite polynomiaL We ob- 
tain, therefore, a non-regular lattice of spins interacting with a strength inversely 
proportional to the square of their distance pH]. The spectrum of this model can 



be found quite easily with the above method. Again, we refer to the literature for 
details [^, |6^ . 



In the continuum limit (A^ -^ oo) the antiferromagnetic version of both the above 
models become c = 1 conformal field theories, the HS containing both chiral sectors 
while the inhomogeneous harmonic one containing just one sector. 

Other models exist and can be solved in this spirit: hierarchical (many-coupling) 
models [0, supersymmetric models ||6^, 'twisted' models |6^ etc. All, however. 



work only for the fundamental representation of some internal group. The big, 
important open problem is to crack a particle system with a higher representation 
for the spins and arbitrary coupling strength. If this is done, through the freezing 
trick we will be able to solve a spin chain with spins in a higher representation. This 
is interesting since we could then see if the antiferromagnetic system for integer 



SU{2) spins develops a mass gap, according to the Haldane conjecture [[70 |. 



6 Exclusion statistics 

So far we approached the problem of statistics in a 'fundamental' way, trying to give 
a reasonable definition and presenting systems that realized this definition. In this 
last section we shall give a 'phenomenological' approach, based on the state-counting 
properties of a system whose dynamics may remain, otherwise, undetermined. It 
is based on a notion already familiar from the GSM model, the notion of 'state 
repulsion' or 'exclusion.' This will lead to Haldane's 'exclusion statistics.' 

6.1 Motivation from the CSM model 

Although Haldane derived his definition from the properties of the HS spin chain, 
we will use the CSM model instead. Consider the 'principal' quantum number of 
this model, that is, the one in terms of which the energy eigenvalues are a sum of 
independent terms: quasiexcitation numbers rij for the Calogero and quasimomenta 
Pi for the Sutherland model. They obey an 'enhanced' exclusion principle where 
nearby values can be no closer that i units. It seems as if each of them occupies 
i places in the single-particle Hilbert space, instead of one. So, let us define the 
dimensionality d{N) of the Hilbert space of states available to an additional particle 
given that there are already A^ particles in the system. (Some high-energy cutoff 
is needed, of course, to make this finite.) The 'exclusion statistics parameter' g is, 
then defined as WT 

Ad , , 

AA^ ^ ^ 

If this parameter is independent of A^, or becomes a constant for high enough A^, 
then we say that the system obeys exclusion statistics of order g. Clearly g = 0,1 
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corresponds to bosons and fermions, respectively. From the discussion of the CSM 
model we would conclude that it obeys exclusion statistics with g = i. 

Note that, in principle, this definition applies to systems of arbitrary dimension- 
ality. The fact that only one- dimensional exclusion statistics systems have been as 
yet identified presumably points to something essentially one- dimensional in this 
definition. Note also that g can be fractional: d{N) could be either ill-defined (as, 
in fact, in the CSM model) so that it can assume "approximate" fractional values, 
or it could depend on A^ in a way that g becomes constant only for large enough 
AA^, like, e.g., d{N) = [N/2] with [] the integer part, giving g = 1/2. 

6.2 Semiclassics — Heuristics 



Before examining the consequences of ( |199|) let us make some heuristic semiclassical 



arguments about phase space volume and exclusion to give more substance to it. 
Please view the following discussion simply as additional motivation -do not take it 
seriously! 

Semiclassically, the number of states of a system is given by the volume of its 
phase space in units of h for each pair of canonical variables. Let us consider for 
simplicity the minimal space {q,p)- It could correspond, for instance, to the relative 
coordinate and momentum of two particles on the line. The standard lagrangian 
would look like 

L = pq- H{p, q) (200) 

If we want the presence of the one particle to 'knock out' g states out of the 
Hilbert space, we should include some term in the classical action that, effectively, 
reduces the volume (in this case, area) of phase space by gh. This area is given as 



A = / dpdq = i pdq (201) 

Jd JdD 

where ^ is a domain in phase space and dD its boundary. If the domain does not 
include the point q = p = the area should be the standard one, since we are talking 
about a region of phase space where the particles are apart. If, on the other hand, 
q = p = is included, the particles are together and the area should diminish by 
gh = 27cg (taking h = 1). 

This reminds us of the Aharonov-Bohm effect. There must be a 'Dirac string' 
piercing the point (0, 0) in phase space giving this extra contribution when circled. 
So we must add to the action a term 



pdq — qdp f ( \q^ 

2 , x2 2 = -9 d[ atan— 
p"^ + X'^q'^ •' \ P 



-^9 I J, ,; 2 =-9 I d[ atan-^ ) (202) 



where A is any positive constant. This amounts to adding to the lagrangian the 

extra term 

pq-qq 
' = -^^ p2 + A2g2 (203) 
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This extra contribution is a total time derivative (a topological term) so it will not 
change the equations of motion. It is expected, though, that it will change the 
quantum mechanical states as described above. 

Let us consider the above in the specific example of two particles in an external 
harmonic oscillator potential. After separating the center of mass phase space, the 
lagrangian for the relative part becomes 

. pq-qq 1 2/ 2 , 2n fonA\ 

Lg=pq-ug-^^^-^--to {p +q) (204) 

For later convenience we chose \ = uj in the g-ieiva. The effect of this term is to 
shift the Poisson brackets between p and q. We will follow the simple approach to 
canonical quantization of defining new 'polar' variables p, 9 as 

p + iuq= v^pe*^ (205) 

In terms of these the lagrangian becomes 

L = C-p'~g)e-]^up' (206) 

so the canonical momentum of 9 is 

T^e = \p'-9 (207) 

and the hamiltonian becomes 

H = uj{iTe + g) (208) 

Quantum mechanically the operator 9 is not quite well defined. In the absence of g, 
it would correspond to the 'phase' of the annihilation operator which cannot be a 
well-defined hermitian operator: e*^ should be unitary, yet it decreases the eigenvalue 
of p^ = a^a by one unit, which is impossible for the ground state. Nevertheless, we 
will proceed qualitatively and see what we get in this case. 

9 being a phase, its momentum vre should be quantized to integer values. Since 
the particles are identical, the change of relative coordinates p — > — j9, q -^ —q should 
also be a symmetry. This corresponds to ^ -^ ^ + vr. Choosing the wavefunction to 
transform as ± itself under this shift, which amounts to quantizing the particles as 
bosons or fermions, further restricts the values of T^g to even or odd integers. Finally, 
since p^ is a positive definite operator, we have 

^p2 = TT, + (7 > (209) 

(clearly the problems with the definition of 9 are hidden in the above constraint). 

Choosing the even eigenvalues for -kq (bosonic case), and for g not greater that 
2, the spectrum of the relative hamiltonian becomes 

E = uj{2n + g) (210) 
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for n a nonnegative integer. This is the excitation spectrum of the Calogero model! 
To see this, add the center of mass oscillator energy Ecm = ^f^t for m a nonnegative 
integer and define 

ni = n, n2 = m + n + g (211) 

Then the full spectrum becomes 

E = uj{ni + 77-2) with ni < n2 — g (212) 

that is, the excitations of the Calogero model in terms of the quasiexcitation numbers 
ni^2 satisfying the 'least distance t constraint ioi C = g. 

This will serve as enough motivation that the idea of phase space exclusion should 
be related to eigenvalue repulsion a la CSM model. The questions about the ground 
state energy, g > 2 etc. that are left hanging are simply set aside -this model is not 
treated seriously here. If anyone is interested, of course, they are welcome to polish 
it! 

6.3 Exclusion statistical mechanics 

We will implement exclusion statistics in terms of the possible quantum numbers 
of A^ particles placed in K single-particle states. Arranging the K states in a lin- 
ear fashion and implementing the 'least distance g^ constraint for the particles, for 
integer g we get []7T], |7^ 

[K-{g-l){N-l)]\ 
^^^' ^^ = Ar![i^-^(Ar-l)-l]! ^^13) 

possible combinations for allowed values for the particle quantum numbers (to prove 
it is a simple combinatorial matter). Extrapolated to arbitrary (fractional) g, this 
defines the state multiplicity of A^ particles placed in K states. Clearly g = 0,1 
reproduces the Bose and Fermi result. 

This can be used to derive the statistical mechanical properties of a grand en- 
semble of such particles. Considering all K states to be at nearby energy e and 
maximizing the Gibbs factor 

D{K, N)e^^''-'^^ = maximum (214) 

in terms of A^ for fixed temperature T = l/kj3 and chemical potential n we obtain 



for the distribution n = N/K for large A^, /sT g, |T|, |73 

n= A^, —^ = ePi^-^) (215) 

1 + (^-1)m; ' (l-w)9 ^ ' 

The expression for n cannot be found analytically except for a few values of g. 

It is obvious from the above that fractional statistics particles are not 'indepen- 
dent', in the sense that one cannot derive their statistical mechanics by considering 
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each single-particle state as an independent system to be filled by particles. (Op- 
erator constructions that incorporate some of the features of exclusion statistics do 
exist 1?^ but they further require specific interactions, else they produce Gentile 
statistics.) There is no microscopic formulation: we need to start from K states and 
take the limit of large K. Further, it is not clear how we could implement exclusion 
statistics for an interacting system, where we cannot use single-particle states as a 
convenient basis. In the following we will sketch how the above difficulties can be 
overcome 



lb 



The starting point will be the grand partition function for particles in K states 
of like energy 

oo 

Z{K,z) = J2 D{K,N)z^ (216) 

N=0 

where z = expP{fi — e) is the fugacity. Z{K,z) should be extensive, which means 
that, for large K it should become the K-th power of a function of z 

1 1 oo 

lim —lnZ{K,x)= lim — In V /^(iT, A^)z^ = In V Pn^" (217) 

In the above we introduced the quantities P„ which play a role analogous to the 
allowed states of occupation of a single level in the cases of fermions and bosons: 
P„ = 1 for bosons, while Pq i = 1, P„>i = for fermions. We can call them 'a priori 



probabilities of occupation' |]T5| or 'fractional dimensionality of states' [^ according 
to taste. 

That Pn exist should be guaranteed from the proper extensive behavior of the 
exclusion statistical mechanics. The problem is to calculate them. For this, we make 
a technical trick: to derive the multiplicity of states D{K, N) above we placed the K 
states on a line. Let's place them on a circle instead, and implement the 'minimum 
distance g^ rule there. That shouldn't infiuence the statistical behavior at large K. 
This modifies the combinatorics into 

We can check that D{K,N) leads to the same statistical distribution n ( pi5| ) as 
D{K,N). We now notice that Z{K,z) defined in terms of D{K,N) becomes an 
exact K-th power for all K. We can, then, calculate Pn from D{l,n). The result is 

Pn=f[(l-'-) (219) 



m=2 



m 



Unless g = 0,1, Pn are fractional and always become negative for some values of n 
||77[ ] . So their interpretation as probabilities or space dimensions must be taken with 
a grain of salt. They are, at any rate, useful tools for describing these systems. 
The single-level grand partition function Z{z) can be shown to satisfy 

Z9 - Z^-i = X (220) 
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from which the corresponding relation (|215D for n = zd^ In Z follows. 



Remember that the CSM system enjoyed a particle-hole coupling-inverting du- 
ality. Since exclusion statistics are directly extracted from the properties of these 
systems, we anticipate a similar duality here [^, ^ . Indeed, we can show that the 
grand partition function satisfies 

^ ^ 1 (221) 



Z{g,x-9) Z{g-\x) 

Note that, interestingly, a similar relation holds for Gentile statistics of maximum 
occupancy p. The grand partition function obviously is 

Z(p, z) = l + z + z\.. + zP = ^ (222) 

z — I 

Identifying the maximum occupancy p of Gentile statistics states with the parameter 
1/(7, the above expression satisfies ( ^21| ) although Z(p, z) is not a priori defined for 
fractional p. 

From the distribution function n{g, z) ( pl5| ) we can derive in the standard way 
the low-temperature Sommerfeld expansion for the energy 



E{T) = I deep{e)n{g, e) = E{0) + ^ Cr^{g)T^+'E^ (223) 

where p(e) is the density of states and /i is fixed through 

N= fdep{e)n{g,e) (224) 



En are T- and (^-independent energy integrals. The coefficients C„ essentially de- 
termine the low-temperature heat capacity of the system. The are calculated as 



73 

Co = (third law of thermodynamics) (225) 

Ci = — (same as in conformal field theory) (226) 

C2 = 2C(3)(1-^) etc. (227) 

where ({x) is the Riemann zeta function. In general, C„ is a polynomial of degree 
n — 1 in g. 

The grand potential G = — kT In Z can be expressed in terms of cluster coeffi- 
cients Wn'- 



ln^P„z" = ^^;2" (228) 



n=0 n=l ^ 



We find for Wr 



m 



n (1 - 5^) (229) 
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m=l 



which is remarkably similar to P„ (except for the range of m). Similar results 



were obtained in |^ for the case of anyons in a strong magnetic field. This is not 
surprising: the lowest Landau level becomes, essentially, a two-dimensional phase 
space, so this is yet another realization of fractional statistics in an effectively one- 
dimensional space. 

Note that the connection between P„ and w„ is the same as the one between 
'disconnected' and 'connected' diagrams in field theory. This, in view also of the 
discussion of the first section, will lead us to a path-integral representation for the 
partition function of exclusion statistics particles in an arbitrary external potential. 

6.4 Exclusion statistics path integral 

We start from the usual euclidean path integral with periodic time (3 for A^ particles 
with action the sum of A^ one-particle actions; we further sum over all particle 
numbers A^ with chemical potential weights e^^'^ /N\ to obtain the grand partition 
function. Since the particles are identical we include a symmetry factor of 1/A^!, but 
we must also sum over paths where particles have exchanged final positions. Thus 
the path integral for each N decomposes into sectors labeled by the elements of the 
permutation group Perm{N). (See the discussion in section 1.) Such permuted 
sectors will be summed with appropriate extra weighting factors, to be determined 
in the sequel. 

By the usual argument, the Gibbs grand potential (the logarithm of the grand 
partition function) will be given by the sum of all connected path integrals. It is 
obvious that these are the ones where the final positions of the particles are a cyclic 
permutation of the original ones, since these are the only elements of S'tv that cannot 
be written as a product of commuting elements. Clearly all diagrams corresponding 
to different cyclic permutations are equal. There are {N — 1)! such permutations. 
So, overall, this diagram will carry a factor of ((A^ — 1)!/A^! = 1/A^. This is nothing 
but the 'symmetry factor,' familiar from Feynman diagrams, corresponding to cyclic 
relabelings of particle coordinates (compare, already, with the factors 1/n included 
in dH)). 

Cyclic permutation diagrams being connected, they only contain one 'thread' of 
particle worldline. They really correspond to one particle wrapping A^ times around 
euclidean time and thus will reproduce the single-particle partition function with 
temperature parameter N(3. So, in the grand partition function they will contribute 
the terms proportional to e'^'^^e^'^^ = z^ (the first factor is the single-particle 
Boltzmann factor for temperature parameter Nj3 and the second is the chemical 
potential factor). 

But these are the terms appearing in (|228D for each energy level! The extra 



weighting factors that must be included for these diagrams in order to reproduce 
exclusion statistics are, then, identified to the cluster coefficients wn- For g = Q^l 
we have wn = 1 and wn = (— )^~^, respectively, which is, indeed, the correct factor 
implied by Xr{P) as discussed in section 1, where R is the trivial (bosons) or the 
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antisymmetric (fermions) irrep of S^ and P a cyclic permutation. 

In conclusion, if we weight these configurations with the extra factors wn we will 
reproduce the grand potential of a distribution of g-ons on the energy levels of the 
one-body problem, that is [|75[ 



fi(/3,/x) = Y: ^'"'^ / e-i:^^i^-[^"(*")] n Dxn{K) (230) 

Af=l ^^ •' n=l 

where Se is the one-particle euclidean action and the paths obey the boundary 
conditions 

Xn{P) = X„+i(0) , x^iP) = Xo(0) (231) 

(x can be in arbitrary dimensions.) 

By exponentiating, the full A^-body partition function will be the path integral 
over all disconnected components, with appropriate symmetry factors and a factor 
of Wn for each connected n-particle component. More directly, we can omit the 
symmetry factors, sum explicitly over all equivalent permutations in each sector 
and divide by A^!: 

oo ^N r N ^ 

^(/^>/^) = E ITTtE W^P / e-S.^i^-[-(*")l n Dx^itn) (232) 

7V=1 -'*• p -^ n=l 

where the paths obey the boundary condition 

{a;i(/5), . . . xn{P)} = P{xi{0), . . . xn{0)} (233) 

The weighting factor Wp for each permutation depends only on the conjugacy class 
of the permutation and is calculated through w„ as 

Wp = Wnr--Wn, (234) 

where ni, . . . ra^ are the cycles of P. 

It is clear that the above path integral is not unitary, since the weights Wn are 
not phases and they do not provide true representations of the permutation group 
(unlike the ^' = 0,1 cases). This does not matter: it is simply used as a tool to 
derive the statistical mechanics of exclusion statistics. We are not going to calculate 
and propagators or other processes with it. 

Once we have the path integral (|232|) we can easily extend the notion of exclusion 
statistics to interacting particles: we simply replace the action J2n SeIxu] by the full 
interacting iV-particle action, thus circumventing all difficulties with combinatorial 
formulae. In the interacting case one has to work with the full grand partition 
function, rather than the grand potential, since topologically disconnected diagrams 
are still dynamically connected through the interactions and do not factorize. 

In conclusion, we have a path-integral way of defining exclusion statistics for 
interacting systems through appropriate symmetry factors included in each permu- 
tation sector (analogous to S{P) of section 1), calculated via ( ^281 ). There is, yet. 



no application of this procedure to an interesting interacting situation. 
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6.5 Is this the only 'exclusion' statistics? 

After giving so much motivation for it, the question seems almost offensive! But let 
us examine it for the sake of ensuring that we have the full picture. 

Exclusion statistics in terms of 'repulsion rules' for the quantum numbers is, 
certainly, as we defined it. But we started from ( |199|) which defines statistics in terms 
of the Hilbert space dimension of an additional particle in the system. Consider, 
then, a chunk of K levels in the CSM model {K quasimomentum values in the 
Sutherland model or K quasiexcitation levels for the harmonic Calogero model). 
Place A^ particles there, and see what space is left for an additional particle. 

If the A^ particles are packed 'densely' (like, e.g., in the Fermi sea-like ground 
state) they certainly take up a piece of approximately IN spaces. So the levels left 
for an extra particle are d = K — iN and we get g = —Ad/ AN = i as expected. 

What if, however, the A^ particles are 'sparse' in K; that is, if the distances 
between them are all bigger than, say, 2£? Then each particle makes unavailable i 
states either way around it, so, overall, the available space for the extra particle has 
diminished by {2i — 1)N. In this situation we would get g = 2£ — 1. Clearly this 
definition implies different statistics for different situations (dense or dilute). 

Can we define another statistics that matches closer the Hilbert space definition? 
Clearly any choice of combinatorial formula for D{K, N) that has the right extensive 
properties is an alternative definition of some statistics. Alternatively, each choice 
of Pn or Wn amounts to some statistics. Let's make the simplest choice: 

Wn = (-«)"-' (235) 

which, in the path integral, corresponds to one factor of —a for each unavoidable 
particle crossing. This leads to the statistical distribution -n 

n = -, — -, (236) 



which was analyzed in [30| as the simplest imaginable generalization of the Fermi and 



Bose distribution. The combinatorial formula for D{K, N) for the above a-statistics 



IS 



^^' _ K{K-a){K-2a)---{K-{N-l)a) 

N\{^ -N)\ ~ m' 



D = «" .ruK .n, = ^ ^TT^ ^ — (237) 



This can be thought as a realization of the exclusion statistics Hilbert space idea: 
the first particle put in the system has K states to choose, the next has K — 
a due to the presence of the previous one an so on, and dividing by A^! avoids 
overcounting. Fermions and bosons correspond to a = 1 and a = —1 respectively, 
a = corresponds to Boltzmann 'classical' statistics. This should also be clear from 
the path integral: no configurations where particles have exchanged positions are 
allowed, since their weighting factor is zero, but factors of 1/A^! are still included. 
The corresponding single-level P„ are 

Pn - n ^ (238) 

mil 1 + "^ 
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For a = 1/p with p integer (a fraction of a fermion), the above 'probabihties' are all 
positive for n up to p and vanish beyond that. For a < all probabilities are positive 
and nonzero. Thus, the above system has a bosonic {a < 0) and a fermionic {a > 0) 
sector, with Boltzmann statistics as the separator. It is a plausible alternative 
definition of exclusion statistics, and has many appealing features, not shared by 
the standard exclusion statistics, such as positive probabilities, a maximum single- 
level occupancy in accordance with the fraction of a fermion that a represents, and 
analytic expressions for all thermodynamic quantities. 

This should demonstrate that the route to alternative definitions, with nice fea- 
tures too, is open. The real issue is whether a dynamical system realizes these other 
statistics, and, at this point, this is not clear. 



7 Epilogue 

The topic of statistics will not cease fascinating at least a few physicists working 
in various fields. It challenges our understanding of the fundamentals of quantum 
mechanics and matter, expands the mathematical tools of the trade, promises new 
results in specific systems and is fun to think about. 

It is anyone's guess whether considerations coming purely out of statistics will 
produce a (major or minor) breakthrough in any physics problem. A safer predic- 
tion, however, is that the idea of generalized statistics will re-emerge in different 
contexts as we strive to understand new physical systems and develop useful frames 
of mind. After all, statistics has barely been touched for objects such as strings 
and membranes, which are increasingly the entities of theoretical choice for the 
fundamental constituents of the universe. The future promises to be fun! 
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